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a)  SCIENTIFIC  RESEARCH  GOALS 

>■  The  main  objective  of  this  research  is  to  investigate  the 

relationships  between  modern  concepts  o-f  nonlinear  dynamical 

systems,  chaos  theory  and  -fully  turbulent  boundary  layers.  The 

approach  that  we  are  taking  is  to  search  -for  a  connection  between 

strange  attractors  and  the  large-scale  coherent  structures  which 

appear  to  play  a  dominant  role  in  the  dynamics  o-f  turbulent  shear 

flows.  ,  Within  the  -framework  of  this  research  contract  we  are 

attacking  the  high  Reynolds  number  turbulent  boundary  layer  from 
♦ 

two  fronts:  lx )  From  the  lower  Reynolds  number  transitional 
boundary  layers,  which  have  the  proper  flow  geometry  but  a  lower 
Reynolds  number  ana  therefore  lower  dimensional  dynamics;  and 
(ii)  from  the  high  Reynolds  number  free  shear  layer  flows,  which 
have  a  simpler  flow  geometry  than  the  fully  turbulent  boundary 
layer,  but  for  which  the  coherent  structures  and  their  dynamics 
aopear  to  be  much  simpler.  ;  Because  of  the  somewhat  reduced 
complexity,  investigating  these  flows  first  will  offer  a  far 
better  opportunity  to  establish  a  connection  of  transitional 
and/or  turbulent  behavior  with  low-dimensional  strange 
attractors.  With  tue  experience  and  understanding  gained  from 
these  simpler  flows,  we  will  then  be  in  a  position  to 


the  fully  turbulent  boundary  law=*^ 


successfully  attack 


7  w* 
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b  >  MAJOR  progress  in  the  past  year 

Continuation  of  build— up  and/or  modification  of  various 
facilities  for  investigating  coherent  structures  and  their 
connection  to  chaos: 

i)  Shear— layer  facility  (water)  with  extensive  use  of 
various  flow  visual ization  techniques  and  advanced 
measurement  techniques  (LDV,  particle  tracking,  etc.). 
(A.  Glezer) 

ii)  Jet  Tacility  (air).  (R.  A.  Petersen) 

iii)  Wind  tunnel  (air)  for  wake  studies.  (F.  Champagne,  I. 
Wygnanski > 

iv)  New  water  channel  for  investigation  of  boundary  layer 
transition.  (H.  Fasal ,  A.  Glezer,  R.  A.  Petersen) 

v)  New  versatile  water  channel  for  wakes  (plane  and 
ax i symmetric) ,  shear  layers,  etc.  (F.  Champagne,  A. 
Glezer,  R.  A.  Petersen,  I.  Wygnanski) 

“:l°®CLGl§yiital_  Investigations 

-  Connection  between  correlation  dimension  and  large  coherent 
structures  in  jet  -'lows  via  the  use  of  diagnostic  tools.  (R. 
Petersen,  I.  Wygnanski,  H.  Fiedler.) 

—  Application  of  the  proper  orthogonal  decomposition  technique  to 
a  harmonically  forced  plane  mixing  layer.  (A.  Glezer  and  A 
Pearl  stein) 

ItlS:2C.e  t  ic  a L  Work  in  Progress 

—  Dynamics  and  Stability  of  Soli  ton  Solutions  of  the  Damped  and 
Driven  Sine-Gordan  Equation.  (D.  McLaughlin,  A.  Pearlstein) 

I 
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-  Development  o-f  Numerical  Methods  for  Determining  Dimensions, 
Attractors  from  "Real"  Data  (data  from  experiments,  or  Navier — 
Stokes  simulations),  (S.  Lichter,  D.  Rand,  A.  Newell) 

-  Coherent  Mode  -  Stochastic  Mode  Decomposition  of  Turbulent 
Fields,  (A.  Newell) 

-  Development  of  Navier — Stokes  Methods  for  Investigating  Dynamics 
of  Coherent  Structures  in  Free  Shear  Layers  and  Wakes.  (H. 
rase! 5 

-  Development  of  Navier — Stokes  Methods  for  Simulating  Later 
Stages  of  Transition  (appearance  of  random  motion).  (H,  Fasel ) 

-  Nonperiodic  Flow  Generated  by  an  Oscillating  Two-Dimensional 
Cylinder.  (A.  Pearlstein) 

-  Development  of  chaos  diagnostic  tools  that  incorporate  spatial 
information  (i.e. ,  spatial  coherence)  and  application  to 
transitional  and  turbulent  jets.  (D.  Rand,  J.  Caputo,  R. 
Petersen) 


c)  PLANS  FOR  NEXT.  YEAR  IS  RESEARCH 

All  the  work  listed  under  b  is  in  progress  and  will  be 
continued  in  the  coming  year. 

Newly  added  activities  will  be: 

-  Chaos  diagnostics  and  analysis  of  numerical  data  obtained  from 


Navier  Stokes  simulation,  of  laminar-turbulent  transition  in 
boundary  layers.  (A.  Slezer,  A.  Pearlstein,  H.  Fasel) 


Devel  opiT.ent: 


numerical  model  for  simulations 


;i  ti  or 


and  turbulence  in  axisymmetric  jets  (to  allow  direct  comparison 


with  jet  experiments) „  (H.  Fasel  and  R.  Li) 


Phase  space  statistics  analysis  of  data  from  jet  experiments 


(R.  Petersen) 
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Cd,e,f,g)  List  of  Publications,  Reports,  Presentations 

List  of  Major  Pub j_icatjgn s /Reports /Presentations  of  Participants 
in  Thfs  Research 

i.  Papers  Published' in  Re-fereed  •lournals: 


H,  Fasel 

"Generation  o-f  Three-Dimensional ,  Spatially  Amplified 
Disturbance  Waves  by  Periodic  Flowing  and  Suction  in  a  Flat 
Plate  Boundary  Laver-"  with  U.  Konzelmanr.  and  l!.  Rist),  ZAMM 

67,  p. " 293,  < 1 937 > . 

"The  Influence  of  Wall  Temperature  on  the  Development  of 
Tol lmien-Schl ichting  Waves  in  Boundary  Layer  Flows,"  (with 
H.  Bestek  and  P.  Dittrich),  ZAMM  67,  p.  256,  (1987). 

"Breakdown  of  a  Two-Dimensional  Laminar  Separation  Bubble  in 
a  Flat  Plate  Boundary  Layer,'*  O'ith  K.  Gruber)  ,  ZAMM  67.  p. 

286,  (19e7). 

"Numerical  Investigation  of  the  Three-Dimensional 
Development  in  Boundary  Layer  Transition,"  (with  LL  Rist  and 
Konzelmann) ,  AIAA  Pager  87- f 203 ,  (1987). 

"Interaction  Between  a  Tollmien-Schlichting  Wave  and  a 
Laminar  Separation  Bubble,"  (with  K.  Gruber  and  H.  Bestek), 

AIAA  Pager  B7-1256,  (1987).  ■ 

"Nonlinear  Interactions  of  Two-Dimensional  Tollmien- 
Schlichting  Waves  in  a  Flat  Plate  Boundary  Layer,:  (with  J, 

Currie  and  A.  Thumm) ,  ZAMM  68  (to  be  published  in  1988). 

"Numerical  Investigation  of  the  Stability  of  Taylor  Vortex 
Flow  in  a  Wide  Gap,"  (with  J.  A.  Meyer),  ZAMM  68  (to  be 
published  in  1988), 

"Numerical  Simulation  of  Subharmonic  Resonance  in  Boundary 
Layer  Transition,"  (with  U.  Konzelmann  and  U.  Rist),  ZAMM 

68,  (to  be  published  in  1988). 

"Numerical  ’  Simulation  of  the  Boundary  Layer  Transition 
Process  Near  the  Secondary  Instability,"  (with  U.  Rist  and 
U.  Konzelmann),  ZAMM  68,  (to  be. published  in  1988). 

"Numerical  Investigation  of  the  Gnset  of  Chaos  in  the  Flow 
Between  Rotating  Cylinders,"  (with  E.  Laurien) ,  ZAMM  68  (to 
be  published  in  1988). 

A-  Glezer 

"On  the  Formation  of  Vortex  Rings,"  accepted  Phvs_>_  Fluids. 

"On  the  Breakdown  of  the  Wave  Packet  Trailing  a  Turbulent  Spot  in 
a  Laminar  Boundary  Layer*"  (with  Y.  Katz v  and  I.  J.  Wygnanski), 
accepted  Fluid  Mechanics. 


"On  the 
Chiarchia 


D.  McLaughlin 
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Weak  Limit  of  Rapidly  Oscillating  Waves,"  (with  L. 
and  N.  Ercolani),  to  appear  in  Duke  Math^  CL  ,  1988. 

"Geometry  o-f  the  Modulational  Instability'.  Part  Is  Local 
Results;  Part  2:  Global  Results,"  (with  N.  Ercolani  and  M.  6. 
Forest-)  submitted  to  Comm^  Pure  Aggl^  Mat  h^ ,  1983. 

Part  1  also  Memoirs  of  the  AMS,  1988  and  Part  2  Memoirs  of  the 
AMS  1988  and  Physica  D  1988  "Connection  Between  Homoclinic 
Structure  of  Isospectral  Set,  in  Instabilities. 

"Geometry  o-f  the  Modulatianal  Instability,  Part  3:  Homoclinic 

Orbits  -For  _he  Period  Sine-Gordon  Equation,"  (with  N.  M.  Ercolani 
and  M.  G.  Forest)  to  appear  Physica  D,  1988. 

"A  Quasi -periodic  Route  to  Chaos  in  a  Near  Integrable  F‘DE," 
(with  A.  Bishop,  M.  G.  Forest,  and  E.  Overman),  Physics  Let t^. 
Vgli  126_.  bBjl  335-340,  1983. 

“Stability  and  Bifurcation  of  Time-periodic  Solutions  of  the 
Damped  and  Driven -Sine-Gordon  Equation,"  (with  A.  Pearlstein  ano 
G.  Terrones)  preprint,  Univ.  of  Aris.,  1988. 

"Stability  and  Bifurcation  of  Spatially  Coherent  Solutions  of  the 
Damped-Driven  NL3  Equation,"  (with  G.  Terrones  and  A.  Pearlstein, 
in  preparation. 

"Notes  on  Melnikov  Integrals  for  Models  of  the  Driven  Pendulum 
Chain,"  (with  N.  M.  Ercolani  and  M.  G.  Forest)  preprint,  Univ.  of 
Aris. ,  1988. 

"Chaos  in  a  Perturbed  Sine-Gordon  Equation  and  in  a  Truncated 
Modal  System,"  (with  A.  R.  Bishop,  R.  Flesch,  M,  G.  Forest,  and 
E.  A.  Overman  II)  preprint,  Los  Alamos,  1987. 

"Modal  Representations  of  Chaotic  Attractors  for  the  Driven 
Damped  Pendulum  Chain,"  (with  A.  R.  Bishop,  M.  G.  Forest,  and  E. 
A.  Overman  II)  preprint,  Ohio  State  University  ,  1988. 

"Finite  Amplitude  Modal  Equations  for  Nearly  Integrable  PDE's," 
(with  N.  M.  Ercolani  and  M»  G.  Forest)  preprint,  Ohio  State 
University,  1938. 

A.  Newell 

"Fixed  Points  and  Chaotic  Dynamics  of  an  Infinite  Dimensional 
Map,"  (with  J.  W.  Moloney,  H.  Mdachihara,  and  D.  W.  McLaughlin, 
Noise  and  Fractals,  pp.  137-186,  1937. 

"A  Calculus  Curriculum  for  the  Nineties,-"  (with  D.  Lovelock) 
ECSE.e'ldi.ngs  NRC-MAA  "Calculus  Curr i.cul.umii  Oct.  i987. 

"Snell's  Laws  at  the  Interface  Between  Nonlinear  Dielectrics, " 
(with  A.  B.  Aceves  and  J.  V.  Moloney)  ,  F-hvsi_cs  Letters  Aj_  Val_ 
129,  No.  4,  1988,  pp.  231-235. 
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"Re-flection  and  Transmission  a-F  S.el-f -focused  Channels  at 
Nonlinear  Dielectric  Interfaces,  (with  A.  B.  Aceves  and  J.  V. 
Maloney) . 

"Tra jectori es  of  Surface  Waves  at  the  Interface  Between  Nonlinear 
Dielectrics  (with  A.  Aceves  and  J.  Moloney),  submitted  to  Physics 
Letters . 

"Wavenumber  Selection  of  Convection  Rolls  in  a  Box , "  (with  W. 
Arter  and  A.  Bernoff ) ,  accepted  Physics  of  Fluids  Letters. 

"Numerical  Simulation  of  Hayleigh-Benard  Convection  in  Shallow 
Tanks,  (with  W.  Arter;  ,  accepted  Physi.cs  of  Fluids. 

"Solitary  waves  as  Fixed  Points  of  Infinite-Dimensional  Maps  for 
an  Optical  Bistable  Ring  Cavity:  Analysis,"  (with  H.  Adachihara, 
D.  W.  McLaughlin,  and  J.  V.  Moloney),  to  appear  J.  of 
L'^themati caf  Physics. 

"Lax  Pairs,  Backlund  Transformations  and  Special  Solutions  for 
Ordinary  Differential  Equations,"  (with  J.  D.  Gibbon,  M.  Tabor, 
and  Y.  Zeng)  ,  to  appear  Ngnl.ineari.ty . 

Turbulent  Transport  and  the  Random  Occurrence  of  Coherent 
Events,  (with  D.  Rand  and  D.  Russell) ,  to  appear  in  special 
volume  of  Phvsica  Dj.  Nonlinear  Phenomena,  Sept.  1938. 

"Theory  of  Beam  Reflection,  Transmission,  Trapping,  and  Breakup 
at  Nonlinear  Optical  Interfaces,  (with  A.  B.  Aceves  and  J.  V. 
Moloney) ,  to  appear  Ogtical  Bistability  I  1938. 

Turbulent  Dissipation  Rates  and  the  Random  Occurrence  of 
Coherent  Events,"  (with  D.  A.  Rand  and  D.  Russell).  to  appear 
Phy^-i  Letters  A . 


A.  Pearl  stein 

Convection  in  Variable  Viscosity  Fluids:  An  Assessment 
of  Approximate  Viscosi ty-Temperature  Relations,"  Phvsics  of 
Fluids^  31,  pp.  1380-1385,  1983. 

"Stability  of  Free  Convection  Flows  of  Variable  Viscosity  Fluids 
in  Vertical  and  Inclined  Slots,"  (with  Y,-M.  Chen)  accepted  J. 
Fi.ui.cJ  Mechanics. 

"Temperature  Distributions  in  a  Laminar-Flow  Tubular 
Photoreactor , "  (with  F.  Chen)  accepted  AIChE  Journal. 

"Low  Peclet  Number  Heat  7r ansf er  in  a  Laminar  Tube  Flow  Subjected 
to  Axially  Varying  Wall  Heat  Flux,"  (with  B.  P.  Dempsev)  accepted 
JL  Heat  Jr ansf er . 

"Efficient  Transf ormation  of  Certain  Singular  Polynomial  Matrix 
Eigenvalue  Problems,"  (with  D.  A.  Goussis)  accepted  J. 
Cgmgutatfgnal.  Physics, 

"Removal  of  Infinite  Eigenvalues  in  the  Generalized  Matrix 
Eigenvalue  Problem,"  (with  D.  A.  Goussis)  accepted  J. 
SSH'SyLa  t  j.  ona  1  Physics .  — 
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2.  Technical  Reports: 


H.  Fasel 


QNR  Contract  N00014— 85— K— 0412  yearly  repcrt,  Mar.  19.87 


3.  Presentations 
a.  Invited 


.  H.  Fasel 

"Navier-Stokes  Simulations  of  Laminar-Turbul ent  Transi  ti  on  ,  " 
International  Symposium  on  Computational  Fluid  Dynamics, 
Sydney,  Australia,  14-17  August  1987. 

"Numerical  Investigation  o-f  Stability  and  Transition  in 
Shear  Flows,"  Euromech  Colloquium  on  Boundary  Layer 
Instability  and  Transition,  Exeter,  England,  21-25  Sept. 


"Numerical  Investigation  o-f  Transition  in  Wall  Bounded  Sheer 
Flows,"  Cali-f.  Inst,  o-f  Technology,  April  1838- 

"Numerical  Simulation  o-f  Breakdown  to  Turbulence,"  Illinois 
Institute  of  Technology,  Chicago,  Ill-,  May  1988. 

"Stability  and  Transition  and  Transition  Control  in  Shear  Flows," 
McDonnell  Douglas  Research  Labs,  St.  Louis,  Mo.,  May  1988. 

A.  Glezer 

DARPA/URI— Brawn/Vale  Conference  on  Turbulent  Structures  in  Free 
Shear  Flows  and  Their  Detection  by  Proper  Orthogonal 
Decomposition,  Newport,  RI,  June  1988. 

Cornell  University,  Mechanical  and  Aerospace  Engineering  Seminar, 
March  1988. 

Princeton  University,  Mechanical  and  Aerospace  Engineering 
Seminar,  Dec.  1987. 

D.  McLaughlin 

Joint  Summer  Research  Conferences  in  the  Mathematical  Sciences, 
University  of  Colorado,  1937. 

Second  Howard  University  Symposium  on  Nonlinear  Semigroups, 
Partial  Differential  Equations,  and  Attractors,  1987. 

Summer  School,  2  weeks  of  lectures,  Ravello,  Italy,  1987. 

Math  Colloquium,  University  of  California,  Davis,  1988. 


University  of  California,  San 


Nonlinear  Science  Colloquium, 

Diego,  1988. 

Nonlinear  Schroedinger  Con-fersnce,  France,  19SS. 

"Critical  Level  Sets  o-f  Integrable  PDE's,"  Solitons  in  Physics 
and  Mathematics,  Minneapolis,  Minn.,  Sept.  1988. 

A.  Pearlstein  . 

Gordon  Research  Conference  on  Oscillations  and  Dynamic 
Instabilities  in  Chemical  Systems,  Plymouth,  NH,  -July  1988. 

DARPA/UR I -Brown /Yale  Conference  on  Turbulence  Structures  in  Free 
Shear  Flaws  and  Their  Detection  by  Proper  Orthogonal 
Decomposition,  Newport,  RI,  June  1988. 


University  of  California,  San  Diego,  Dept,  of  Applied  Mechanics 
and  Engineering  Science,  Feb.  1988. 

Yale  University,  Applied  Mechanics  Colloquium,  Nov.  1987. 
b.  "  Contributed 


H.  Fasel 

"Nonlinear  Interactions  of  Two-Dimensional  Tollmien- 
Schlichting  Waves  in  a  Flat  Plate  Boundary  Layer,"  (with  J. 
Currie  and  A.  Thumm) ,  SAMM  Conference,  Stuttoart,  Aoril 
1987. ■ 


"Numerical  Simulation  of  Subharmonic  Resonanc  e  in  Boundary 

Layer  Transition,"  (with  U.  Konzelmann  and  'J.  Rist),  GAMM 
Conference,  Stuttgart,  April  1987. 

“Numerical  Investigation  of  the  Onset  of  Chaos  in  the  Flow 
Between  Rotating  Cylinders,"  (with  E.  Laurien) ,  GAMM 
Conference,  Stuttgart,  April  1987. 

"Numerical  Investigation  of  the  Stability  of  Taylor  Vorte:< 
Flow  in  a  Wide  Gap,"  (with  J.  A.  Meyer),  GAMM  Conference, 
Stuttgart,  April  1987. 


" Numerical  Simulation  of  the  Boundary 
Process  Near  the  Secondary  Instability," 
U.  Konzelmann) ,  GAMM  Conference,  Stuttgart 


Layer  Transition 
(with  U.  Rist  and 
,  April  1987. 


"Interaction  Between  a  Tollmien— Schlichting  Wave  and  a 
Laminar  Separation  Bubble,"  (with  K.  Gruber  and  H.  Bestek) , 
AIAA  19th  Fluid  Dynamics,  Plasma  Dynamics  and  Lasers 
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DETAILED  DISCUSSION  OF  SOME  RESEARCH  IN  PROGRESS 
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Th«  available  data  consists  of  10  chaneia  of  U(x,y,t) 
Before  Applying  the  P.O.D. 

*  Smoothing  in  time 

*  Piecewise  polynomial  fit  in  Y 
a  Subtract  the  mean  velocity 
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APPLICATION  OF  THE  P.O.D.  TO  THE  DATA 
Lum I  ay (1967) ;  Aubry,,  Holmes,  Lumley  &.  Stone  (1986) 

Previous  applications  have  been  to  fully  developed  turbulent 
flows  that  are  : 

*  Statistica I ly  Homogeneous  in  Two  Directions 

*  Statistca I ly  Stationary 


DESCRIPTION  OF  THE  P.O.D. 


U(i,i)  . ; . u(i,n) 

A  *  . ;o  <>;  45 ; 

uCm»i) . !!!!!!!!!!!!!!!  !u(m|n) 


Singular  Value  Decomposition  of  A  '(SVD)  is  equivalent  :;o  compu- 

T 

ting  the  eigen  values  of  A«A 

The  eigen  values  measure  energy  or  action  in  correponding  eigen 
vectors  (  spatial  modes  )  . 

Eigen  vectors  constitute  an  orthogonal  basis  for  representing 
U(*,y,t)  . 
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Conclusions 

•  Can  apply  POD  "to  nonsPaAionar^ 

nonhomo^enous  PurbulenP  -flows 

•  Can  repcessnp  -flovo  by  supeeposiV-con 
o-C  small  number  of  POD  modes 
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Can  OSS  small  pqrp  of  "  temporal 
record  Po  qe.oera.-Ve  POO  modcss.- 
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Turbulent  Dissipation  Rates 
and  the  Random  Occurrence  of  Coherent  Events 

by 

Alan  C.  Newell,  David  A.  Rand*,  David  Russell  f 
Department  of  Mathematics 
University  of  Arizona 
Tucson,  AZ  85721 

Abstract 

In  this  letter,  we  suggest  that  the  transport  properties  and  dissipation  rates  of  a 
wide  class  of  turbulent  flows  are  determined  by  the  random  occurrence  of  coherent 
events  that  correspond  to  certain  orbits  which  we  call  homoclinic  excursions  in  the 
high  dimensional  strange  attractor.  Homoclinic  excursions  are  trajectories  in  the 
non-compact  phase  space  that  are  attracted  to  special  orbits  which  connect  saddle 
points  in  the  finite  region  of  phase  space  to  infinity  and  represent  coherent  structures 
in  the  flow  field. 


Introduction 

It  is  generally  accepted  that  turbulence  in  shear  flows  at  Reynolds  numbers  of 
103  and  higher  has  a  large  number  of  active  degrees  of  freedom.  Although  estimates 
of  Lyapunt  and  Hausdorff  dimensions  ax'*  strictly  upper  bounds,  the  rapid  loss  of 
spatial  correlations  over  distances  of  the  Tu,  lor  microscale  R~?  and  the  broadband 
wavenumber  spectrum  suggest  that  many  modes  are  playing  an  active  role  in  the 
dynamics.  The  R*  estimate  [1],  which  is  consistent  with  the  intuitive  idea  of  Landau 
that  it  is  necessary  to  resolve  a  turbulent  flow  field  in  a  box  of  volume  V  down  to 
the  Kolmogoroff  inner  scale  of  (r/3/s)*  (1/  is  kinematic  viscosity,  e~v  fv  |  V  v.\2dx 
is  the  energy  dissipation  rate  and  is  independent  of  u  for  a  large  range  of  large 
Reynolds  numbers),  means  that  if  R  is  104,  R*  is  a  billion!  It  is  unlikely  that 
replacing  the  Navier-Stokes  equations  by  a  system  of  a  billion  o.d.e.s.  will  either 
bring  much  insight  into  the  nature  of  turbulence  or  make  it  possible  to  calculate 
the  invariant  measure  on  the  attractor  needed  to  compute  averages.  Therefore, 
whereas  the  concept  of  a  strange  attractor  on  which  the  flow  is  everywhere  unstable 
and  ergodic  is  a  valuable  and  necessary  one,  a  new  idea  is  needed  if  one  is  to  be 
able  to  compute  in  a  practical  way  average  flow  quantities  and  in  particular  those 
averages  which  represent  transport  properties,  the  flux  of  heat  across  a  convection 
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f  Present  Address:  Center  for  Nonlinear  Studies,  Los  Alamos,  NM  87545 


1 


layer,  momentum  across  a  boundary  layer,  angular  momentum  between  cylinders 
rotating  at  different  velocities,  mass  flux  down  a  pipe  or  the  amount  of  heavier  fluid 
which  falls  through  a  lighter  one.  It  is  also  important  to  calculate  the  dissipation 
rate  which,  in  systems  where  energy  is  fed  in  at  low  wavenumbers  and  removed 
at  high  ones,  is  equivalent  to  computing  the  flux  or  cascade  of  energy  across  any 
middle  wavenumber.  One  of  the  main  goals  of  a  turbulence  theory  is  to  predict  the 
transport  and  the  dissipation  rate  as  a  function  of  the  applied  stress.  The  purpose  of 
this  paper  is  to  suggest  one  new  avenue  of  approach  to  these  problems  of  transport 
in  physical  and  wave  space. 

The  basic  idea  is  fairly  simple.  There  is  some  evidence  that,  in  a  variety  of 
situations,  transport  is  associated  with  organized  flow  structures.  If  this  is  the 
case,  one  would  like  to  identify  and  isolate  those  parts  of  the  attractor  and  those 
orbits  in  phase  space  which  are  the  principal  contributors  to  a  particular  flux.  In 
this  letter  we  focus  on  a  class  of  such  special  orbits  which  connect  saddle  points  in 
the  finite  part  of  phase  space  to  infinity  and  which  are  the  principal  contributors 
to  the  dissipation  rate.  The  ideas  are  first  illustrated  in  terms  of  a  simple,  but 
nontrivial,  example  ,  the  two  dimensional  forced,  damped,  nonlinear  Schrodinger 
(NLS)  equation 

i/>,-iV*il>-iW\2rl,  =  F-D  (1) 

which  arises  in  models  of  Langmuir  turbulence.  Although  this  model  is  a  poor 
approximation  to  the  full  Zakharov  equations  on  several  counts,  it  nevertheless 
retains  the  essential  feature,  emphasized  for  many  years  by  Soviet  colleagues  [2], 
that,  in  the  limit  of  strong  ion  damping,  dissipation  is  mainly  due  to  the  collapse 
of  filaments  rather  than  energy  transfer  to  high  wavenumbers  by  resonant  wave- 
wave  interactions.  The  filaments  are  closely  related  to  an  exact,  singular,  localized 
solution  of  the  unperturbed  and  conservative  NLS  equation  (1),  whose  modulus  has 
the  shape 

\ip(x, t)|  =  A ^R(rj)  ,  V  =  A-1|z  —  xQ\  (2) 


near  the  blow-up  point  z0.  In  (2),  R(rj)  is  the  unique  solution  of  R"  +  — 

R  +  R3  —  0,  R'( 0)  =  0,  i?(oo)  =  0  without  zeros  in  (0,  oo)  and  A (t)  is  a  time 
dependent  function  which  becomes  zero  after  a  finite  time  to*  The  rate  of  collapse 
A(i)  has  not  been  yet  satisfactorily  determined.  The  two  most  recent  efforts  [3] 
both  agree  that  it  has  be  form  /(t)(t0  —  where  /(<)  tends  to  zero  very  slowly 
(e.g.(ln(i0  —  i))-1  or  (lnln(t0  —  i))-1)  but  differ  on  the  exact  asymptotic  answer. 
A  sufficient  but  not  necessary  condition  for  collapse  is  that  the  Hamiltonian  H  = 
j*(i  VVM2  —  \ \ib\^)dx  is  negative.  The  negativity  of  3  guarantees  that  =  f  |V»|2dz, 
also  a  constant  of  the  motion  for  the  NLS  equation,  is  greater  than  the  minimum 
power  p  =  2ir  J0°°  R2rjdTj  =  0.29  needed  to  sustain  collapse.  If  =  p,  Weinstein 
[4]  has  shown  that,  if  the  solution  blows  up  in  finite  time,  then  rp(x,  t )  converges 
to  (2)  in  the  H 1  norm.  On  the  other  hand,  if  <  p,  global  existence  is  assured 
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and  no  blow-up  occurs.  Our  interest  is  in  larger  values  of  n$.  In  particular,  an 
initial  condition  which  is  a  small  perturbation  of  the  unstable,  monochromatic  state 
ip  =  exp (ik  •  x  +  i(k 2  —  |2)t)  in  a  square  box  of  side  length  L  with  k  = 

2vL~l(m,  n),  m,  n  integers  (m  —  n  =  0  corresponds  to  the  x  independent  state)  for 
which  —  \ipo\2L2  will  develop  approximately  n^p-"1  collapsing  filaments,  each 
having  precisely  the  shape  (2)  near  its  collapse  site  xo  and  time  to*  In  the  phase 
space  of  the  unperturbed  system,  each  collapsing  filament  will  be  represented  by 
an  orbit  which  joins  the  saddle  point  representing  the  unstable  solution  to  infinity. 
We  give  these  orbits  a  special  name,  heteroclinic  connections  to  infinity  or  HCI’s, 
because  all  orbits  which  escape  to  infinity  are  attracted  to  them  even  though  the 
unperturbed  system  is  conservative.  La  other  words,  all  solutions  of  NLS  which 
become  singular  in  finite  time  have  the  shape  (2)  as  they  approach  infinity.  Of 
course,  in  the  damped  system,  the  singularity  is  never  quite  reached.  Once  the 
filament  diameter  is  of  the  order  of  the  Debye  radius,  Landau  damping  (represented 
by  I?  in  (1))  becomes  important  and  the  energy  is  transferred  from  the  wavefield 
to  electrons,  a  process  know  as  burnout.  The  process  of  dissipation  changes  the 
attracting  point  at  infinity  to  a  repeller  and  the  system  returns  to  the  finite  part 
of  phase  space  where  it  again  comes  under  the  influence  of  a  saddle  point  and  the 
cycle  is  repeated.  We  call  these  large  excursions,  from  the  finite  part  of  phase  space 
to  the  neighborhood  of  infinit-'”  and  back  again,  homoclinic  excursions.  They  are 
organized  by  the  HCI’s  in  the  sense  that  as  the  trajectory  approaches  infinity,  the 
corresponding  solution  in  physical  space  approaches  a  very  special  shape  and  in 
particular  the  amount  of  energy  it  carries  off  to  the  dissipation  cemetery  at  infinity 
is  known.  Before  attempting  to  mathematize  these  notions  any  further,  we  present 
some  numerical  evidence  to  support  our  picture. 


Numerical  Results 

We  now  turn  to  the  main  results  of  the  numerical  study.  We  integrated  (1)  using 
a  split-step  Former  algorithm  on  a  grid  of  (128)2  points.  The  time  step  used  was 
10~4  and  adequate  to  resolve  the  largest  linear  (k2)  and  nonlinear  (|*/’|2)  frequencies 
encountered  in  our  simulations.  Aliasing  errors  were  judged  to  be  insignificant  by 
making  spot  comparisons  with  a  dealiased  code  on  a  (256)2  grid.  More  details 
are  given  in  [5].  Energy  is  fed  into  the  system  at  low  k  and  removed  at  high  k 
through  F  and  D  whose  Fourier  transforms  are  7o(k)V’fc  and  'yt(k)rpk  respectively 
(  ipk  is  the  Fourier  transform  of  ip(x ,  t)).  The  support  of  7 t,(k)  is  near  k  =  0  and 
the  support  of  7 e(k)  is  beyond  kd  where  kd  is  the  Debye  wavenumber.  It  turned 
out  that  in  most  of  our  simulations  the  global  energy  H  (not  a  constant  of  motion 
for  (I))  was  positive  during  all  time  intervals  in  which  there  were  collapse  events. 
Figure  1  shows  snapshots  of  |^>2(x,t)|  before,  during  and  after  the  collapse.  Figure 
2  displays  the  instantaneous  rate  of  dissipation  7 (i)  =  2  X)£7e(&)|1/,fc|2,  its  integral 
T(<)  =  —  f  and  the  global  maximum  of  \ip\2  over  a  time  interval  including 

several  collapses  in  a  strong-turbulent  regime  after  transients  have  died  out  and 
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the  system  has  reached  a  statistically  steady  state.  — r  is  just  the  energy  lost  due 
to  dissipation  as  a  function  of  time.  Observe-  that  a  large  fraction  of  the  energy 
dissipation  occurs  in  sudden  jumps  and  is  directly  correlated  to  the  collapse  events. 
Energy  n ^  —  J^k  iV’fcl2  is  considered  to  be  dissipated  due  to  collapse  only  if  7  >  70, 
where  70  is  a  threshold  parameter.  If  this  inequality  holds  over  the  interval  [ii,  *2]} 
then  the  energy  burnt  out  is  T(ti)  —  Tfa),  and  the  total  energy  lost  due  to  collapse 
for  the  entire  simulation  is  the  sum  of  all  such  differences.  Denote  this  total  energy 
lost  due  to  collapse  by  Snc  and  the  total  energy  dissipated  by  Sn.  We  have  plotted 
the  ratio  Snc/6n  as  a  function  of  70  for  several  different  turbulent  regimes.  As  70 
decreases  the  ratio  increases  and  is  of  course  one  at  70  —  0.  70  =  1  was  the  smallest 
value  that  gave  ratios  judged  to  be  free  from  contamination  by  non-collapse  events 
for  all  cases  studied,  so  it  was  chosen  as  the  standard  cut-off  rate.  With  this  choice 
of  70,  as  much  as  80%  of  the  total  energy  dissipated  is  dissipated  by  collapsing 
filaments  in  the  most  energetic  case  examined  (<  >=  10).  In  the  weakest 

case  (<  riff/  >=  5)  more  than  70%  of  the  dissipated  energy  was  lost  in  burned  out 
filaments.  Furthermore,  these  estimates  of  the  energy  dissipation  are  conservative 
because  we  have  not  yet  included  the  more  gradual  loss  due  to  the  decay  of  the 
high  k  remnant  left  over  from  the  collapse.  The  reason  for  this  remnant,  peculiar 
to  the  nonlinear  Schrodinger  equation,  is  that  a  singular  filament  of  this  equation 
carries  the  minimum  threshold  energy  p  required  to  sustain  collapse.  (In  contrast, 
the  energy  of  a  collapsing  Zakharov  filament,  again  a  constant  of  the  motion,  can 
take  on  a  continuous  range  of  values  above  the  critical  threshold.)  Therefore,  when 
a  significant  portion  of  the  filament’s  spectral  energy  lies  in  the  dissipative  range 
(k  >  kd),  some  of  its  energy  is  lost,  collapse  is  arrested  and  only  a  partial  burnout 
occurs.  It  leaves  behind  a  remnant  in  the  form  of  broadened  concentric  cylindrical 
shells  of  field  energy  centered  at  the  collapse  site  (see  Figure  3).  In  order  to  follow 
the  energy  associated  with  the  remnant,  we  monitored  the  field  energy  inside  a  small 
cell  centered  at  the  collapse  site  and  observed  that  the  burnout  of  the  central  part 
of  the  filament,  or  core,  causes  the  fastest  depletion  of  energy  in  the  cell,  but  also 
that  it  is  immediately  followed  by  a  slower  depletion.  The  slower  loss  is  primarily 
due  to  the  fluxing  of  the  remnant  through  the  cell  boundary.  If  we  add  the  energy 
loss  rp,0  <  r  <  1,  due  to  the  burnout  of  the  core  to  that  carried  out  of  the  cell 
by  the  remnant,  we  obtain  a  total  loss  of  energy  in  the  cell  very  close  to  the  total 
amount  y  carried  by  the  collapsing  filament.  Figure  4  shows  the  disribution  of  r. 

The  remnant  plays  two  roles  in  the  dissipation  process.  First,  it  provides  a 
nucleating  center  for  new  collapses,  a  fact  we  have  verified  by  examining  the  spatial 
distribution  of  collapse  sites.  As  a  consequence  the  frequency  of  collapses  increases. 
We  found  that  for  7 j(&)  =  2 6(k  —  1)  (the  beam  drives  just  the  lowest  wavenumber 
in  the  box),  the  distribution  of  collapse  times  is  nearly  Poisson  with  a  mean  time 
between  events  of  w”1  =<  r  >~  .08.  Second,  the  gradual  damping  of  the  remnant 
enhances  the  ambient  dissipation  rate  <  74  >  (i.e.,  the  dissipation  rate  averaged 
over  time  intervals  free  from  collapse)  over  that  amount  we  would  expect  from 
other  dissipative  mechanisms.  In  order  to  verify  this,  we  carried  out  the  following 
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relaxation  experiment.  By  changing  the  sign  of  the  nonlinearity  in  equation  (1)  from 
plus  to  minus,  we  obtain  a  system  with  the  same  linear  dispersion  as  the  plus  system 
and  therefore  presumably  the  same  resonant  wave  interaction  character  but  without 
the  unstable  modulational  and  nucleational  instabilities  (which  will  be  discussed  in 
the  next  section)  and  collapsing  solutions.  Initializing  the  two  systems  identically, 
we  ran  them  with  the  forcing  turned  off  (7 &  =  0)  and  compared  their  dissipation 
rates  between  collapses  of  the  (+)  system,  before  their  energies  had  separated  by 
more  than  10%  of  the  common  initial  value.  The  ambient  dissipation  rate  of  the  (+) 
system  was  typically  50%  greater  than  that  of  the  non-collapsing  (-)  system  while 
the  collapse  remnants  were  dissipated.  Including  collapse  events  in  the  tally,  the 
(+)  system  lost  energy  three  times  as  fast  as  the  (-)  system.  We  conclude,  therefore, 
that  the  amount  of  energy  dissipated  by  non-collapse  events  is  at  most  one-quarter 
of  the  total  energy  dissipated.  Therefore,  the  energy  dissipation  budget  is  as  follows. 
The  average  total  dissipation  rate  <  7  >=  wc  <  r  >  p-f-  <  74  >=  ucp  <  r  > 

since  <  7  >  depends  only  on  the  mean  turbulence  level  <n$  >  but, 

as  we  have  pointed  out,  <  7^4  >  and  wc  depend  both  on  <  >  and  the  nature  of 

the  dissipation.  The  average  function  <  r  >  lost  in  burnout  depends  only  on  the 
latter.  We  plotted  both  <  7  >  and  ue  as  a  function  of  <  >  and  found  that, 

as  expected  from  above,  the  product  <  7  >  is  almost  independent  of  <  >. 

Futhermore  it  was  equal  to  .13  which  is  equal  to  the  product  of  p  (which  was  .29), 
<  r  >  (which  was  .3)  and  (which  was  about  1.3).  Graphs  showing  the  loss  of 
energy  from  a  small  cell  surrounding  the  collapsing  filament,  jjf-  as  a  function  of 
70,  <  7  >  and  wc  as  a  function  of  <  >,  the  distribution  of  collapse  times,  the 

angle-averaged  correlation  function  and  actual  trajectories  illustrating  homoclmic 
excursions  will  be  given  elsewhere  [5]. 


Whereas  these  results  confirm  the  thesis  that  coherent  collapse  events  dominate 
dissipation,  the  expense  involved  in  two  dimensional  experiments  did  not  allow  us 
to  run  a  sufficient  number  of  simulations  in  order  to  determine  in  what  limit,  if 
any,  approaches  unity.  In  connection  with  this  question,  and  in  parallel  with 
the  ideas  of  DiPema  and  Majda  [6],  one  would  also  like  to  develop  the  notion  of 
a  weak  solution  in  which  solutions  of  the  undamped  equations  could  be  continued 
beyond  the  collapse  time  by  simply  deleting  the  collapsed  filament  and  lowering  the 
L 2  norm  of  the  solution  by  a  fixed  amount.  In  order  to  examine  these  questions, 
we  simulated  the  one  dimensional  nonlinear  Schrodinger  and  Zakharov  models  with 
quintic  nonlinearities  on  grids  of  1024  and  256  points  respectively.  Aliasing  errors 
were  removed  by  smoothly  interpolating  ip  onto  a  grid  of  8  x  1024  or  256  points  before 
forming  the  nonlinear  frequency  |^>|4.  The  results  for  NLS  showed  that,  for  the  same 


turbulence  levels,  the  percentage  of  the  dissipatic 


3 tinted  for  by  collapses 


rose  from  59%  at  kd  —  128  (where,  because  the  dissipation  effects  axe  clearly  felt 


in  the  early  stages,  there  are  many  failed  attempts  to  form  collapsing  filaments),  to 
72%  for  kd  =  256,  to  82%  for  kd  =  512.  The  average  loss  of  energy  per  collapse 
decreases  with  kd  (the  amount  of  energy  greater  than  kd  in  the  Fourier  transform 
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of  the  collapsing  filament  decreases  with  increasing  kd )  but  the  frequency  of  events 
increases  proportionately  so  that  the  average  dissipation  rate  remains  the  same. 
The  evidence  clearly  suggests  that  in  the  large  kd  limit,  all  the  energy  dissipated  is 
dissipated  by  collapse  events  which  occur  infinitely  often  with  infinitisemal  losses  of 
energy  per  event. 

A  much  more  striking  result  was  obtained  when  we  ran  the  very  same  experiment 
on  the  Zakharov  model 

ipt  —  i'ff’xx  +  ipifr  —  F  —  eD  (3) 

ptt+2vo  pt-  pss  =(M4)k  (4) 


where  F  and  D  are  as  before  and  2v  o  pt  is  a  convolution  integral  modelling  ion 
damping.  The  ion  acoustic  field  p(x,i)  is  no  longer  slaved  to  the  electric  field 
intensity  (for  NLS,  p  in  the  equation  (3)  is  replaced  by  —  |t/>|4).  Therefore,  during 
collapse,  in  which  the  fields  take  on  a  self-similar  form  close  to  (2)  but  in  which 
the  inertial  acceleration  ptt  is  also  important,  the  cavity  formed  by  the  ion  field 
encourages  total  burnout  of  the  filament  (see  Figure  5).  We  found  that  when  e  — 
1,  kd  =  32,  <  Tty  >=  1.625,  the  average  energy  <  r  >  p  lost  per  event  was  1.95, 
the  average  time  between  events  was  .83  and  was  92%.  (For  the  Zakharov 
model,  the  amount  of  energy  carried  in  the  filament,  a  constant  of  the  motion  for 
the  unperturbed  equations,  can  take  on  a  range  of  values  greater  than  the  threshold 
value  p  =  .43  in  this  case.  Thus  r  >  1).  The  distribution  of  r  rose  sharply  after 
r  =  1,  had  a  mean  of  1.95/.43  and  had  a  relatively  long  taiL  Also,  we  observed  that 
all  the  collapses  occurred  in  9.  nonoverlapping  sites,  were  driven  by  the  nucleational 
instability  (see  the  section  after  next  for  discussion),  and  the  sites  drifted  about  the 
box.  When  we  increased  kd  to  64,  <  7ty  >=  1.8,  the  frequency  of  events  increased 
to  (.6P)_1  (there  were  287  events  in  198  time  units)  and  the  distribution  of  r  came 
closer  to  one.  The  average  energy  lost  per  event  was  1.57  units  and  ^  was  95%. 
(In  contrast,  at  the  same  parameter  values  and  over  100  time  units,  the  NLS  model 
had  972  events  with  an  average  energy  loss  of  .09  and  was  72%.)  Further, 
when  we  decreased  e  to  L  (again  kd  =  64),  the  distribution  of  r  came  closer  to 


threshold  and  was  98%.  Also,  for  e  =  1,  when  kd  was  increased  to  128,  ^ 
was  98%.  This  leads  us  to  conjecture  that  as  kd  increases  and  e  decreases,  the 
distribution  of  r  will  cluster  about  one,  and  that  each  collapsing  event  will  bum  off 
exactly  the  threshold  energy  p  =  .43.  In  this  asymptotic  limit,  the  weak  solution 
of  the  unforced,  undamped  equation  in  which  the  turbulence  eventually  decays,  is 
found  by  simply  removing  the  collapsing  filament  from  the  field  and  reducing  the 
Li  norm  by  p.  (The  weak  limit,  i.e.  y(r,  0)  is  the  weak  limit  of  rp(x,  A)  as  A  — ►  0  if 
f  0  )dx  =  lim  /  <p(x)rp(x,  A  )dx  where  (p  is  smooth,  of  a  collapsing  filament 

is  zero  because  it  oscillates  very  fast  and  its  width  decays  at  a  faster  rate  than  its 
amplitude  increases).  We  will  report  more  details  elsewhere.  We  are  also  currently 
testing  some  ideas  concerning  the  estimation  of  the  frequency  of  events. 
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Coherent  Events 


Our  picture,  then,  is  that  in  this  and  a  wide  variety  of  similar  situations,  the 
dissipation  rate  is  governed  by  the  random  occurrence  of  coherent  events.  In  order 
to  compute  a  good  approximation  to  the  dissipation  rate,  one  need  only  know  the 
following  information;  (i)  the  nature  of  the  hcmoclinic  excursion  (equivalently  the 
spatial  and  temporal  shape  of  the  coherent  structure)  and  the  energy  lost  in  the 
cycle  and  (ii)  the  average  frequency  of  occurrence.  Although  finding  approximations 
to  these  quantities  will  not  in  general  be  an  easy  task,  it  should  be  considerably 
simpler  than  finding  the  complete  statistics  on  the  turbulent  flow.  Further,  this 
picture  suggest  that  the  most  natural  decomposition  of  the  turbulent  field  is  one  in 
which  the  field  is  divided  into  low  dimensional,  organized  structures  (the  collapsing 
filaments)  and  a  high  dimensional  component  consisting  of  chaotic  fluctuations  in 
the  field  during  intervals  when  no  collapse  events  occur.  In  phase  space,  the  high 
dimensional  component  will  cause  the  homoclinic  excursion  to  fluctuate  chaotically 
about  the  HCI  but,  because  in  physical  space  these  fluctuations  represent  chaotic 
behavior  away  from  the  collapse  site,  it  is  reasonable  to  expect  that  their  effects 
can  be  captured  using  low  order  statistics  and  that  they  will  not  appreciably  affect 
the  dissipation  rate. 

One  would  like  to  understand  better  the  low  dimensionality  and  attractive  nature 
of  the  organized,  singular  structures.  Why  is  it  that  the  Only  structures  through 
which  the  rp(x,  t)  field  becomes  singular  are  relatively  simple?  We  have  no  mathe¬ 
matically  rigorous  or  even  plausibly  compelling  argument  but  suggest  that  this  is 
a  very  important  question  to  answer.  Our  own  attempt  at  an  answer  is  incom¬ 
plete  but  stresses  the  importance  of  the  scaling  symmetry  which  the  NLS  (and  Za¬ 
kharov)  equation  enjoys:  if  ip(x,  t)  is  a  solution,  so  is  <p(x,  T )  where  =  A~xt/>,  X  — 
A “xx,T  =  X~2t  which  suggests  the  family  ip  =  A"xi2(A“x|xj)exp(— i\~2t)  of  exact 
solutions.  If  ip  is  to  approach  infinity  as  A-1(i),  then  in  order-  that  the  dispersion 
and  time  derivative  terms  can  balance  cubic  nonlinearity,  x  must  scale  as  A-1  and 
t  as  A-2.  Thus  a  singular  solution  in  which  dispersion  can  balance  nonlinearity 
(spatially  independent  solutions,  some  of  which  can  be  singular,  are  unstable)  de¬ 
mands  scale  symmetry.  If  one  inserts  ip  =  A“x<p(X  =  A-Xx,  T  =  f  \~2di)  into 
the  NLS  as  the  leading  term  in  an  asymptotic  expansion  and  if  one  makes  the  fur¬ 
ther  ansatz  that  <p  =  f(x)expid,  where  9  ~  J  X~2dt+  lower  order  terms,  then  the 
demand  that  the  coefficient  of  A-3  be  zero  gives  a  nonlinear  eigenvalue  problem 
for  /(x)  for  which  there  is  a  unique  isotropic  ground  state  solution,  namely  R(rj). 
However,  the  separation  ansatz  is  not  forced  at  this  level.  One  might  conjecture 
that  if  one  attempts  to  construct  a  uniformly  valid  expansion  for  t/>,  then  both  this 
ansatz  and  the  behavior  of  A(t)  might  result  from  solvability  conditions.  This  is  the 
point  of  view  taken  by  Zakharov  [3].  Despite  the  absence  of  a  rigorous  argument, 
it  would  appear,  nevertheless,  that  symmetry  constraints  have  much  to  do  with  the 
low  dimensionality  of  the  singular  structures. 
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Geometry  of  Phase  Space 


The  dominance  of  filaments  in  the  dissipation  process  suggests  that  the  attractor 
A  contains  two  saddle-like  unstable  sets  M  and  S  and  heteroclinic  orbits  between 
them.  M  contains  saddle  points  representing  modulationally  or  nucleationally  un¬ 
stable  (we  discuss  these  shortly)  solutions  of  (1),  (3)  and  (4).  The  stable  manifold 
of  S  at  infinity  is  low  dimensional  and  corresponds  to  the  idealized  collapsing  fil¬ 
aments  of  the  unperturbed  equation.  The  unstable  manifold  of  M  intersects  the 
stable  manifold  of  S.  The  unstable  manifold  of  S  (the  burnout  due  to  Landau 
damping)  intersects  the  stable  manifold  of  M.  The  speed  of  attraction  and  repul¬ 
sion  at  S  are  governed  by  two  entirely  different  processes  and  axe  not  related.  The 
former  depends  on  the  faster  than  exponential  collapse  while  the  latter  depends 
on  dissipation  and,  for  the  Zakharov  equation,  the  slow  relaxation  of  the  ion  hole. 
Based  on  this  picture  A  can  be  roughly  subdivided  into  two  sets  Ac  and  Ah-  When 
the  system  is  in  Ah  ,  it  is  dominated  by  what  we  all  the  hash  modes  consisting  of  the 
background  field  and  radiation  modes  left  over  from  the  formation  and  collapse  of 
the  coherent  structures.  This  set  will  generally  be  large  dimensional  but  contributes 
little  to  the  dissipation  rate.  On  the  other  hand,  when  the  system  is  in  Ac,  it  is 
dominated  by  the  coherent  structures,  and  in  particular  as  the  saddle  point  S  at  in¬ 
finity  is  approached,  the  trajectories  asymptotically  tend  to  the  special  orbit  which 
we  have  earlier  called  the  HCI.  In  other  words,  the  part  of  the  turbulent  solution 
which  provides  the  principal  contribution  to  the  dissipation  rate  is  low  dimensional. 
After  the  energy  is  removed  by  Landau  damping  (or  in  general  by  whatever  dissi¬ 
pation  process  is  relevant)  the  phase  point  returns  to  the  turbulent  soup  Ah  where 
it  again  comes  under  the  influence  of  the  saddle  points  in  M  and  the  cycle  repeats. 

We  now  turn  to  a  discussion  of  the  structure  of  the  phase  space,  the  nature 
of  the  saddle  point  M  and  its  unstable  manifold.  Whereas  the  perturbations  of 
forcing  and  damping  modify  the  phase  space  structure,  it  is  the  topology  of  the 
phase  space  of  the  unperturbed  system  which  sets  the  stage  for  the  large  homo¬ 
clinic  excursions.  T\  .  properties  of  the  unperturbed  system,  which  is  Hamiltonian, 
are  crucial.  One  is  the  fact  that  its  phase  space  contains  many  fixed  points  and 
periodic  orbits  of  saddle  type  and  separatrices  joining  these.  The  second  crucial 
property  is  the  non-compactness  of  the  constant  energy  surfaces  and  the  existence 
of  heteroclinic  connections  (HCI’s)  which  join  certain  of  the  saddle  points  to  in¬ 
finity.  Whereas  Hamiltonian  perturbations  of  this  .system  can  destroy  certain  tori 
through  resonances  and  create  new  chains  of  elliptic  centers  and  hyperbolic  saddles 
and  non-Hamiltonian  perturbations  can  turn  the  centers  into  attracting  sinks  or 
destroy  them  altogether,  the  important  point  to  make  is  that  the  hyperbolic  nature 
of  the  original  saddles  and  the  attracting  nature  of  their  HCI’s  remain  intact  and 
play  a  dominant  role  in  the  dynamics  of  the  perturbed  system.  The  instabilities 
which  give  rise  to  collapsing  filaments  in  (1)  are  closely  related  to  the  naturally 
occurring  instabilities  of  the  unperturbed  NLS  equation,  an  infinite  dimensional, 
Hamiltonian  system  with  Hamiltonian  E  and  the  additional  motion  constants 
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and  J  =  i  J(ip  V  V’*  V’*  V  There  are  two  types  of  unstable  critical  points. 

The  first  are  nonlocal  solutions  which  are  constant  or  purely  periodic  in  space  with 
a  time  behavior  depending  on  amplitude  (so  they  are  really  unstable  periodic  or¬ 
bits).  These  states  axe  unstable  to  sideband  disturbances  (fluctuations  with  a  spatial 
period  close  to  that  of  the  unstable  solutions),  an  instability  known  as  the  modu • 
lational  instability.  The  second  type  axe  called  cavitons  and  axe  analogous  to  the 
solitons  of  the  one  dimensional  equation.  Like  the  collapsing  filament,  they  have 
zero  H  value  and  any  perturbation  which  decreases  S  will  render  them  unstable. 
Because  of  the  existence  of  unstable  solutions,  the  phase  space  of  the  unperturbed 
Hamiltonian  system  contains  many  saddle  points.  Increasing  the  size  of  the  do¬ 
main  or  the  amplitude  of  the  ij>  field  (the  important  parameter  is  the  L2  norm)  will 
create  saddle-center  bifurcations  and  increase  the  number  of  saddles.  Hamiltonian 
perturbations  such  as  the  addition  of  a  term  V(x)|^|2  to  S  can  create  additional 
island  chains  of  hyperbolic  saddles  and  elliptic  centers  and  in  particular,  if  V(x) 
is  a  localized  potential  well,  these  new  critical  points  can  correspond  to  the  just 
mentioned  localized  caviton  states  [7].  The  presence  of  a  small  amount  of  dissi¬ 
pation  will  destroy  centers,  turning  some  into  sinks,  but  the  saddles  will  persist. 
Other  non-Hamiltonian  perturbations,  such  as  forcing,  can  destroy  the  stable  sinks 
and  destabilize  the  weakly  stable  cavitons  which  axe  localized  in  potential  wells.  In 
practice,  this  scenario  often  arises  [8],  particularly  for  the  Zakharov  equations  and 
in  some  instances  for  the  nonlinear  Schrodinger  equation.  If  ion  damping  is  large, 
the  hole  in  the  ion  acoustic  density  field  remains  after  the  burnout  of  the  electric 
field.  The  hole  acts  as  a  slowly  relaxing  potential  which  can  serve  to  focus  electric 
field  energy  into  a  metastable  caviton  [7].  When  the  hole  amplitude  is  sufficiently 
small,  the  caviton  destabilizes  and  the  field  intensity  and  ion  acoustic  field  again 
form  collapsing  filaments  (see  Figure  5).  We  call  this  the  nucleatioiial  instability. 
The  phenomenon  was  discovered  in  numerical  simulations  by  Doolen,  Russell,  Rose 
and  DuBois  [8]  and  has  also  been  confirmed  in  experiments  by  Cheung  and  Wong 
[9].  It'  is  from  a  combination  of  modulational  and  nucleational  instabilities  that  the 
collapsing  filaments  are  bom.  The  cycle  time  of  the  homoclinic  excursion  in  each 
case  depends  on  (a)  the  time  taken  for  the  phase  point  to  come  under  the  influence 
of  a  saddle  instability  and  (b)  the  growth  rate  of  the  instability.  In  the  NLS  case  it 
is  difficult  to  estimate  the  former.  In  the  Zakharov  case,  however,  this  time,  during 
which  the  phase  point  returns  from  S  back  to  Ah,  can  be  estimated  by  calculating 
how  long  it  takes  for  the  hole  to  relax  [10]. 


Simple  Example  of  a  HCI 

The  existence  of  HCI’s  depends  on  two  crucial  properties  of  the  underlying  un¬ 
perturbed  (undamped,  unforced)  equations;  the  presence  of  saddle  points  in  the 
finite  part  of  phase  space  and  the  non-compactness  of  the  energy  surface  on  which 
they  lie  (or  the  intersection  of  the  level  surfaces  of  the  constants  of  the  motion  of 
there  is  more  than  one).  The  existence  of  saddle  points  and  the  noncompactness 
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of  the  energy  surface  means  that  there  are  many  unstable  solutions  whose  unstable 
manifolds  are  unbounded.  We  do  not  yet  understand  why,  in  infinite  dimensional 
systems,  the  HCI  is  attracting  and  low  dimensional  although  we  can  illustrate  the 
former  property  with  a  simple  finite  dimensional  example.  Consider  the  perturbed 
Puffing’s  oscillator  x"+x  —  2x3  —  F(t)  —  (xx~1)nx-,  where  F(t)  is  small  and  the  ad¬ 
ditional  factor  (xXq1  )n  has  been  added  to  the  damping  term  in  order  to  arrange  for 
it  to  turn  on  only  when  x  is  large.  The  unperturbed,  exactly  integrable  Hamiltonian 
system  represents  a  particle  in  a  qu?rtic  potential  V(x)  =  -(x2  —  x4)  with  critical 
points  x'  =  x  —  0  (a  center)  and  x‘  =  0,x  —1  (saddle  points).  The  unstable 
and  stable  manifolds  emanating  from  x  =1  are  the  non-compact  energy  level 

surfaces  E  —  \  —  |x-2  +  |x2  —  |x4  which  join  the  saddles  to  infinity  through  the 
HCI’s  x‘  =  (x2  —  |)  for  x‘  >  0,x  >  and  x’  =  — (x2—  |)  for  x‘  <  0,  x  <  — The 

other  branches,  the  stable  manifolds  of  x‘  and  x  —t.  are  repellers  for  increasing 

time.  All  energy  level  surfaces  E  >  \  and  those  for  E  <  \  for  which  x2  >  -  escape 
to  infinity  in  finite  time  t0  and  it  is  easy  io  show  that  all  escaping  orbits  are  asymp¬ 
totically  equal  at  infinity.  By  this  we  mean  that  the  distance  between  all  escaping 
orbits  tends  to  zero  as  they  approach  infinity.  To  see  this  easily,  consider  the  Lau¬ 
rent  series  about  the  finite  blow-up  time  to  in  the  unperturbed  Duffing  equation; 

x(t)  s=  (t— to)-1  +6(t— 10)+  ’2  c2j+\(E)(t  —  to)2j+1.  Observe  that  the  information 


as  to  which  escaping  orbit  one  is  on  is  contained  in  the  terms  cubic  and  higher  in 
(t— 10).  The  singular  part  of  all  escaping  orbits  are  the  same.  One  can  also  introduce 
a  canonical  change  of  coordinates,  u  —  — x_1 ,  v  —  x2(x'  —  x2  +  |),  in  which  one  can 
see  that  the  orbit  x'  =  x2  —  x  >  — corresponding  to  E  =  §■,  attracts  all  other 
orbits  escaping  to  infinity.  The  fact  that  Hamiltonian  systems  can  have  attracting 
orbits  (the  point  at  infinity  itself  is  not  even  a  fixed  point)  at  infinity  does  not  seem 
to  be  generally  known  although  it  is  apparent  in  the  work  of  Bogoyavlenskii  [11] 
who  examines  the  neighborhood  of  infinity  for  several  Hamiltonian  systems. 


Nonlinear  Schrodinger  HCPs 

This  simple  example  illustrates  how  the  topology  of  the  phase  space  of  the  un¬ 
perturbed  system  is  important  in  controlling  the  dynamics  of  the  perturbed  sys¬ 
tem.  Likewise  in  the  weakly  perturbed  nonlinear  Schrodinger  equation,  in  which 
forcing  is  applied  at  small  wavenumbers  and  the  damping  at  large,  it  is  the  topol¬ 
ogy  of  the  phase  space  structure  of  the  unperturbed  problem  which  dominates  the 
dynamics  of  the  perturbed  problem.  The  unperturbed  problem  is  again  Hamil¬ 
tonian  although^  in  this  case,  it  is  infinite  dimensional  with  motion  constants 
n$  =  f  ij)ij)mdx,  J  =  i  J  r}>  V  V’"  —  V  $)dx  and  E  =  ~  /(vV*  *  “  |t^2^*2)dx. 

In  particular,  the  intersection  of  the  level  surfaces  and  in  particular  the  energy  is 
non  compact  and  there  are  orbits  (HCI’s)  in  the  surface  >  nc  =  p  which  join  the 
saddle  points  representing  modulational  and  nucleational  instabilities  to  infinity. 
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Our  picture  of  the  phase  space  then  is  as  follows.  For  low  levels  of  forcing  (either 
F  is  small  or  the  domain  area  is  small),  the  modulus  of  the  electric  field  grows 
till  >  nc  at  which  stage  the  phase  point  eventually  comes  close  to  an  unstable 
saddle  point  which  is  joined  to  infinity  via  a  collapsing  filament.  The  instability  sets 
in,  the  filament  is  formed,  it  collapses,  a  certain  portion  of  its  energy  is  dissipated 
(depending  on  the  damping  structure)  and  a  remnant  of  high  k  waves  is  left  over. 
For  the  Zakharov  system,  we  expect  a  more  complete  burnout.  The  system  returns 
close  to  its  original  state  and  the  process  is  approximately  repeated  with  the  initial¬ 
ization  of  the  collapsing  filament  being  due  to  either  modulatioual  or  nucleational 
instabilities.  In  particular,  either  the  remnant  or  the  ion  acoustic  field  can  provide 
a  cavity  for  the  nucleation  of  metastable  cavitons  which  collapse  once  the  support 
disperses  away.  For  larger  values  of  the  applied  stress  (either  larger  forcing  or  larger 
boxes),  the  value  of  the  mean  turbulence  level  is  much,  larger  and  can  be  many 
times  nc.  In  this  case,  spatial  correlations  decay  rapidly  and  many  collapsing  fila¬ 
ments  can  occur  at  different  spatial  locations  although,  since  each  event  is  so  rapid, 
they  will  rarely  occur  in  the  same  time  intervals.  For  these  large  turbulence  levels, 
the  saddle  point  M  represents  a  collection  of  saddle  points  in  the  turbulent  soup 
part  of  the  attractor  Ah,  some  of  which  correspond  to  analogues  of  the  modula- 
tional  instability  of  more  complicated  periodic  shapes  and  others  of  which  are  best 
understood  as  nucleational  instabilities. 


Other  turbulent  situations 


We  mention  here  that  a  similar  situation  obtains  in  the  Euler  equations.  For 
low  values  of  the  applied  stress,  laminar  states  can  be  destabilized  by  identifiable 
instability  mechanisms,  centrifugal  instabilities,  mean  flow  profiles  with  inflexional 
points  and  so  on.  These  instabilities  do  not  go  away  when  the  fluid  is  more  highly 
stressed  and  becomes  fully  turbulent.  Indeed  saddle  points  representing  local  inflex¬ 
ional  and  centrifugal  instabilities  remain  very  much  part  of  the  strange  attractor  of 
the  high  Reynolds  number  Navier  Stokes  equations  and  play  a  large  role  in  trans¬ 
ferring  energy  to  high  wavenumbers  where  dissipation  acts.  In  highly  turbulent 
flows,  the  inflexional  profile  of  the  mean  flow  is  not  sustained  for  all  time  uniformly 
in  space,  but  if  it  is  sustained  long  enough  in  a  local  region,  then  rapidly  growing 
packets  of  three  dimensional  inflexional  instabilities  can  erupt  and  carry  energy  off 
to  the  dissipation  cemetaries.  In  addition,  the  level  surfaces  in  the  phase  space  are 
noncompact.  Moreover,  it  is  not  unreasonable  to  argue  that  the  dissipation  rate 
of  shear  flow  turbulence  arises  not  from  the  wave-wave  interaction  familiar  from 
Fourier  space  cumulant  descriptions,  but  instead  to  suggest  it  is  dominated  by  the 
formation  and  destruction  of  thin  vortex  sheets  or  surfaces  containing  vortex  tubes 
(the  state  S )  in  which  almost  all  of  the  vorticity  is  concentrated.  One  might  conjec¬ 
ture  that  the  initial  formation  of  surfaces  of  vorticity  concentration  follows  singular 
solutions  (either  finite  or  infinite  time)  of  the  Euler  equations  driven  principally  by 
inflexional  instabilities.  Once  formed  (the  system  is  close  to  5),  these  sheets  would 
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be  notoriously  unstable,  to  Kelvin-Helmholz  instabilities  of  the  tangential  velocity 
discontinuities  and  to  Taylor-Gortler  centrifugal  instabilities  of  the  helical  flows  in¬ 
duced  by  vortex  lines  embedded  in  the  curved  sheets.  These  secondary  instabilities 
quickly  transfer  the  energy  to  the  viscous  cemetaries  where  energy  is  dissipated  and 
the  system  is  returned  close  to  a  state  where  approximately  the  same  cycle  can 
repeat. 

Finally,  we  suggest  that  these  ideas  are  not  restricted  to  the  approximation  of 
dissipation  rates.  Other  transport  properties  such  as  the  momentum  exchange  from 
the  plate  to  the  outer  flow  in  a  turbulent  boundary  layer  also  seem  to  be  dominated 
by  “organized  events,”  namely  the  burst-sweep  cycle.  It  is  interesting  that  Aubry, 
Holmes,  Lumley  and  Stone  [12]  have  found  a  type  of  homoclinic  excursion  in  their 
analysis  of  the  “long  wave”  structure  of  a  turbulent  boundary  layer  suggesting  the 
beginning  of  a  burst  sweep  cycle,  although  their  excursion  does  not  lead  to  any 
kind  of  singular  behavior.  One  might  conjecture  that  the  latter  could  result  if 
one  were  to  add  a  short  scale  inflexional  wave  component  which  can  result  from 
such  a  wavepacket  being  phase  locked  to  and  growing  on  the  local  distortion  of  the 
turbulent  mean  profile.  Further,  we  suggest  it  may  be  also  valuable  to  apply  these 
ideas  to  heat  transport  at  high  Rayleigh  numbers  (where  thermal  plumes  world 
play  the  role  of  singular  struct  uresyand  Rayleigh-Taylor  instabilities. 

Acknowledgements 

The  authors  are  grateful  for  a  series  of  extremely  useful  conversations  with  Harvey 
Rose,  Don  DuBois,  David  McLaughlin,  Volodja  Zakharov  and  Sacha  Rubenchik. 
This  work  was  supported  by  ONR  Engineering  grant  N0001485K0412  and  AFOSR 
grant  FQ8671-8601551.  David  Russell  was  partially  supported  by  an  ONR  grant  at 


References 


[1.]  D.  Ruelle,  Commun.  Math.  'Phys.  87,  287-302  (1982).  Commun.  Math.  Phys. 
93,  285-300  (1984). 

C.  Foias,  0.  P.  Manley,  R.  Teman  and  Y.  M.  Treve,  Physi'ca  9D,  157-188  (1983). 

C.  Doering,  J.  D.  Gibbon,  D.  D.  Holm  and  B.  Nicolaenko,  to  be  published  in 
Nonlinearity  (1987). 

J.  M.  Hyman  and  B.  Nicolaenko,  Los  Alamos  Report  LA -UR-85-1556. 

[2.]  A  good  review  article  is  “ Collapse  versus  Cavitons ”  by  A.  M. 
Rubenchik,  R.  Z.  Sagdeev  and  V.  E.  Zakharov.  Comments  on  Plasma 
Physics  and  Control  Fusion  9,  183-206  (1985) 

M.  V.  Goldman,  Reviews  of  Modern  Physics  56,  709  (1984). 

[3.]  B.  LeMesurier,  G.  Papanicolaou,  P.  Sulem  and  C.  Sulem,  Private  Communication 
(1988). 

V.  E.  Zakharov,  Private  Communication  (1988). 

[4.]  M.  J.  Weinstein.  Comm.  Partial  Differential  Equations  11,  545-565  (1986). 

[5.]‘  A.C.  Newell,  D.A.  Rand,  D.  Russell  to  appear  Physica  D  Special  Volume’in  honor 
of  Joseph  Ford  (1988). 

[6.]  A.  Majda,  R.  DiPema,  Private  Communication. 

[7.]  H.  A.  Rose  and  M.  I.  Weinstein,  Physica  30D,  207-218  (1988). 

[8.]  G.  D.  Doolen,  D.  DuBois,  H.  A.  Rose,  Phys.  Rev.  Lett.  54,  804  (1985). 

D.  Russell,  D.  Dubois,  H.  A.  Rose.  Phys.  Rev.  Lett.  56,  838  (1986). 

D.  Russell,  D.  DuBois,  H.  A.  Rose.  Phys.  Rev.  Lett.  60,  581  (1988). 

[9.]  A.  Y.  Wong,  P.  Y.  Cheung,  Phys.  Rev.  Lett.  52  12 22  (1984). 

P.  Y.  Cheung,  A.  Y.  Wong,  Phys.  Rev.  Lett.  55  1880  (1985). 

[10.]  H.  A.  Rose,  Private  Communication. 

[11.]  0. 1.  Bogoyavlenskii,  Methods  of  the  Qua.litc.Hve  Theory  of  Dynamical  Systems  in 
Astro-physics  and  Gas  Dynamics ,  Springer- Verlag,  Berlin  (1981). 

[12.]  N.  Aubry,  P.  Holmes,  J.  Lumley,  E.  Stone,  Cornell  University  Preprint  (1988). 


13 


Figure  Captions 


Figure  1.  Snapshots  of  at  times  (a)  beginning  of  (by  during  and  (c)  after  collapse. 

Figure  2.  A  record  over  a  short  period  of  time  of  the  instantaneous  rate  of  dissipation 
7 (i)  =  2S£7e(fc)|V’(£)!2,  its  integral  3T(f)  =  —  f*  y(t')dt',  and  the  global  spatial 
maximum  of  \ip(k)\2  as  functions  of  time  in  the  strong  turbulent  regime  with 
76  =  2 6(k  -  1). 

Figure  3.  Concentric  shells  which  are  the  remnants  of  two  simultaneous  collapses.  Note: 

simultaneous  collapses  are  rare  and  these  were  initiated  by  a  symmetric  initial 
state. 

Figure  4.  The  distribution  of  r  for  two-dimensional  NLS. 

Figure  5.  The  (a)  nucleation,  (b)  collapse  and  (c)  burnout  of  a  filament  of  the  Zakharov 
equation.  The  overshoot  of  the  cavity  during  burnout  encourages  tetal  dissipation 
of  the  filament. 
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Abstract 

We  suggest  that  the  transport  properties  and  dissipation  rates  of  a  wide  class 
of  turbulent  flows  are  determined  by  the  random  occurrence  of  coherent  events 
which  correspond  to  certain  orbits  which  we  call  homoclinic  excursions  in  the  high 
dimensional  strange  attractor.  Homoclinic  excursions  are  trajectories  in  the  non¬ 
compact  phase  space  that  are  attracted  to  special  orbits  which  connect  saddle  points 
in  the  finite  region  of  phase  space  to  infinity  and  represent  coherer*  structures  in  the 
flow  field.  This  picture  also  suggests  that  one  can  compute  fluxe^  using  a  relatively 
low  dimensional  description  of  the  flow.  A  method  for  extracting  the  organized 
structures  from  a  time-series  is  given  and  provides  a  local  analogue  of  the  notion  of 
Lyapunov  exponents. 
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Introduction 

It  is  generally  accepted  that  turbulence  in  shear  flows  at  Reynolds  numbers  of  103 
and  higher  has  a  large  number  of  active  degrees  of  freedom.  Although  estimates  of 
Lyapunov  and  Hausdorff  dimensions  [1]  are  strictly  upper  bounds,  the  rapid  loss  of 
spatial  correlations  over  distances  of  the  Taylor  microscale  R~  *  and  the  broadband 
wavenumber  spectrum  suggest  that  many  modes  are  playing  an  active  role  in  the 
dynamics.  The  R t  estimate,  which  is  consistent  with  the  intuitive  idea  of  Landau 
that  it  is  necessary  to  resolve  a  turbulent  flow  field  in  a  box  of  volume  V  down  to 
the  Kolmogoroff  inner  scale  of  (u^/e)*  (u  is  kinematic  viscosity,  e  —  v  fv  |  V  u\2dx 
is  the  energy  dissipation  rate  and  is  independent  of  v  for  a  large  range  of  large 
Reynolds  numbers),  means  that  if  R  is  104,  R*  is  a  billion!  It  is  unlikely  that 
replacing  the  Navier-Stokes  equations  by  a  system  of  a  billion  o.d.e.s.  will  either 
bring  much  insight  into  the  nature  of  turbulence  or  make  it  possible  to  calculate 
the  invariant  measure  on  the  attractor  needed  to  compute  averages.  Therefore, 
whereas  the  concept  of  a  strange  attractor  on  which  the  flow  is  everywhere  unstable 
and  ergodic  is  a  valuable  and  necessary  one,  a  new  idea  is  needed  if  one  is  to  be 
able  to  compute  in  a  practical  way  average  flow  quantities  and  in  particular  those 
averages  which  rep.esent  transport  properties,  the  flux  of  heat  across  a  convection 
layer,  momentum  across  a  boundary  layer,  angular  momentum  between  cylinders 
rotating  at  different  velocities,  mass  flux  down  a  pipe  or  the  amount  of  heavier  fluid 
which  falls  through  a  lighter  one.  It  is  also  important  to  calculate  the  dissipation 
rate  which,  in  systems  where  energy  is  fed  in  at  low  wavenumbers  and  removed 
at  high  ones,  is  equivalent  to  computing  the  flux  or  cascade  of  energy  across  any 
middle  wavenumber.  One  of  the  main  goals  of  a  turbulence  theory  is  to  predict  the 
transport  and  the  dissipation  rate  as  a  function  of  the  applied  stress.  The  purpose  of 
this  paper  is  to  suggest  one  new  avenue  of  approach  to  these  problems  of  transport 
in  physical  and  wave  space. 

The  basic  idea  is  fairly  simple.  The\e  is  some  evidence  that,  in  a  variety  of 
situations,  transport  is  associated  with  .organized  flow  structures.  If  this  is  the  case, 
one  would  like  to  identify  and  isolate  those  parts  of  the  attractor  and  those  orbits 
in  phase  space  which  are  the  principal  contributors  to  a  particular  flux.  Our  claim 
is  that  they  are  homoclinic  excursions  analogous  to  homoclinic  orbits  which  occur 
at  random  intervals  on  the  attractor  and  connect  that  large  dimensional  part  of  the 
attractor,  consisting  of  states  of  the  system  which  contribute  little  to  the  flux,  to 
itself.  The  homoclinic  excursions  are  arranged  by  orbits  connecting  the  main  part 
of  the  attractor  to  an  organizing  structure  which  is  represented  in  phase  space  by 
a  generalized  saddle.  Of  particular  interest  to  us  is  the  case  where  this  structure 
is  at  infinity  and  corresponds  to  a  singular  solution.  We  give  the  orbits  connecting 
saddles  in  the  attractor  to  the  singular  structure  at  infinity  the  name  heteroclinic 
connections  to  infinity  (HCI’s).  There  are  also  orbits  from  the  singular  structure 
back  to  the  main  part  of  the  attractor.  The  HCI’s  are  important  because  we  shall 
show  that  there  is  a  sense  in  which  they  are  attracting  even  when  the  system  is 
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conservative.  The  existence  of  HCI’s  depends  upon  two  crucial  properties  of  the 
unperturbed  (undamped,  unforced)  equations:  the  presence  of  saddle  points  in  the 
phase  space  of  the  unperturbed  system  and  the  non-compactness  of  the  energy 
surface  on  which  they  lie  (or  the  intersection  of  the  level  surfaces  of  the  constants  of 
motion  if  there  is  more  than  one).  The  existence  of  saddle  points  in  the  phase  space 
of  the  unperturbed  system  and  the  non-compactness  of  the  energy  surface  means 
that  there  are  many  unstable  solutions  whose  unstable  manifolds  are  unbounded. 
The  HCI,  which  lives  in  the  unstable  manifold,  corresponds  in  physical  space  to  a 
coherent  structure  which  is  usually  localized  in  space,  has  a  finite  lifetime  and  is 
connected  with  a  family  of  exact,  albeit  singular,  solutions  of  the  field  equations. 
At  infinity,  where  a  large  portion  of  the  energy  of  these  solutions  is  in  small  scales, 
the  structure  becomes  unstable  to  dissipative  processes  which  drain  its  energy  and 
relax  the  system  back  to  the  neighborhood  of  its  initial  state,  thus  completing  the 
cycle.  Towards  the  end  of  this  paper,  we  will  discuss  a  prescription  for  identifying 
and  finding  these  structures  using  time  series  data. 

Our  picture,  then,  is  that  in  a  wide  variety  of  situations,  turbulent  transport  is 
achieved  and  the  dissipation  rate  is  governed  by  the  random  occurrence  of  coherent 
events.  In  order  to  compute  a  good  approximation  of  the  dissipation  rate  or  flux 
in  question,  one  need  only  know  the  following  information,  (i)  the  nature  of  the 
homoclinic  excursion  (equivalently  the  spatial  and  temporal  shape  of  the  coherent 
structure)  and  the  energy  lost  in  this  cycle,  and  (ii)  the  average  frequency  of  occur¬ 
rence.  Although  finding  approximations  to  these  quantities  will  not  in  general  be 
an  easy  task,  it  should  be  considerably  simpler  than  accurately  following  the  full 
dynamics.  Obtaining  these  estimates  by  largely  deterministic  means  is  not  intended 
to  disregard  the  fact  that  turbulence  is  a  stochastic  process  in  which  statistical  fluc¬ 
tuations  are  important.  What  it  does  suggest,  however,  is  that  the  most  natural 
decomposition  of  the  turbulent  field  is  one  in  which  the  field  is  divided  into  orga¬ 
nized  motions  represented  by  dominant  orbits  in  the  phase  space  and  statistical 
fluctuations  about  these  orbits.  The  method,  discussed  later,  by  which  we  iden¬ 
tify  the  organized  structures  and  their  dynamics  is  a  variation  of  Karhunen-Loeve 
expansions  (later  used  by  Lumley  [2]  in  the  context  of  turbulence)  which  extract 
from  statistical  data  dominant  shapes  and  forms.  Our  method  is  novel  in  that  it 
is  a  weighted  local  version  of  the  Lumley  method  and  is  therefore  able  to  pick  out 
selectively  relatively  rare  coherent  events  which  long  time  averages  would  lose.  As 
we  have  mentioned,  we  believe  that  in  many  cases  these  dominant  shapes  will  be 
closely  related  to  exact,  singular  solutions  of  the  governing  equations  which  act  as 
asymptotic  attractors  in  a  sense  to  be  described  below.  Further,  with  this  decompo¬ 
sition,  it  may  even  be  possible  to  approximate  the  full  dynamics  on  the  attractor,  at 
least  to  some  degree  of  accuracy,  by  a  finite  and  low  dimensional  system  of  o.d.e.’s, 
with  stochastic  coefficients  reflecting  the  influence  of  the  many,  many  other  active 
degrees  of  freedom  contained  in  the  fluctuations  about  the  principal  orbits.  It  is  our 
thesis  that  a  good  approximation  to  the  transport  may  be  obtained  either  by  ig¬ 
noring  altogether  or  by  simply  taking  average  values  of  these  stochastic  coefficients. 
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In  any  event,  we  suggest  that  using  this  decomposition,  the  fluctuations  can  be  de¬ 
scribed  with  low  order  statistics.  The  principal  difference  with  previous  thinking  [3] 
in  which  one  attempts  to  embed  the  attractor  in  a  large  but  finite  dimensional  Eu¬ 
clidean  manifold  (the  inertial  manifold  approach)  is  that  we  no  longer  demand  that 
the  approximation  to  the  turbulent  field  keeps  the  phase  point,  which  represents 
the  state  of  the  system,  always  in  the  attractor.  Rather  we  simply  demand  that  it 
moves  parallel  to  the  real  attractor,  respecting  the  directions  associated  with  large 
homoclinic  excursions  which  we  claim  are  largely  responsible  for  transport.  This  re¬ 
quirement,  much  weaker  than  the  cone  condition  of  inertial  manifold  theory,  means 
that  far  fewer  coordinates  are  needed  to  describe  those  particular  features  of  the 
dynamics  associated  with  the  homoclinic  excursions. 

Before  attempting  to  mathematize  these  notions  any  further,  let  us  first  consider 
a  representative  case  study  in  which  these  ideas  take  concrete  form. 

The  Zakharov  model  of  Langmuir  turbulence. 

The  Zakharov  equations, 

3  oj 

V2(*Vt  +  i7e°^+  -uprj  V2v>)  \P(V9  +  E0)  +  V  1  -F,  (1) 

Ptt  +  2vi  opt-c23S72p=  — v2  I  V  <f>\2  (2) 

couple  the  density  fluctuations  p(x,  t )  and  the  longitudinal  component  of  the  electric 
field  E  =  |  V(^eWf’<+(*))+Eo  in  a  plasma  and  constitute  a  popular  nonlinear  fluid 
model  of  Langmuir  turbulence  [4].  Here  up  is  the  plasma  frequency  at  which  a  cold 
plasma  is  resonant  in  the  longitudinal  mode  and  the  envelope  (f>(x ,  t)  is  assumed  to 
vary  slowly  in  time  compared  to  w”1.  The  linear,  dissipative  convolution  opera¬ 
tors  7e  and  7,  are  used  to  model  Landau  damping;  both  increase  with  increasing 
wavenumber  but  7e  does  so  abruptly  in  order  to  cut  off  the  coherent  longitudinal 
oscillations  as  the  Debye  wavenumber  k<i  —  —  is  approached;  r<f  is  the  distance  at 
which  the  Coulomb  potential  of  a  bare  charge  is  diminished  exponentially  by  the 
ambient  plasma.  F  is  a  forcing  term  that  injects  electric  field  energy  into  the  plasma 
and  is  also  modelled  by  a  linear  convolution  operator  7 j  0  E.  Equation  (2)  is  the  ion 
acoustic  wave  equation  in  which  ion  density  fluctuations  are  driven  by  gradients  of 
the  electric  field  intensity,  the  so-called  pondermotive  force.  If  the  ions  are  quasi¬ 
stationary,  the  dominant  term  on  the  left  hand  side  of  (2)  is  the  last  and,  in  this 
case,  the  ion  density  fluctuation  p  is  directly  slaved  to  the  field  intensity  |  V  ^*l2-  In 
one  spatial  dimension  and  with  suitable  rescalings,  the  Zakharov  equations  reduce 
to  the  familiar  nonlinear  Schrodinger  equation  in  dimensionless  variables  for  the 
electric  field  envelope  ?/>(£,  t) 

^  +  iy2il>  +  i\ip\2il>  =  F-D.  (3) 


4 


Our  case  study  will  principally  focus  on  the  two  dimensional  version  of  (3).  Although 
this  model  is  a  poor  approximation  to  the  real  physics  on  several  counts  (it  ignores 
finite  sound  speed  effects  and  the  role  that  sound  waves  play  in  energy  transfer,  it 
assumes  the  ion  density  is  slaved  to  the  electric  field  intensity,  it  is  two  dimensional 
and  the  premises  from  which  it  is  derived  become  suspect  as  the  electric  field  reaches 
its  singular  limit),  nevertheless  it  retains  the  essential  features  of  the  full  dynamics 
in  the  sense  that  the  dissipation  rate  is  controlled  by  the  filament  dynamics  Energy 
is  fed  in  at  low  k  and  removed  at  high  k  through  F  and  D  whose  respective.  Fourier 
transforms  axe  7 b(k)ip(k)  and  7e(k)’0(&)(i/>(£,  t)  is  the  Fourier  transform  of  ip(x,t)) 
with  the  support  of  7 i(k)  being  near  k  =  0  and  the  support  of  7 e(k)  being  k  >  kd. 

For  many  years  it  has  been  conjectured  by  Soviet  colleagues  [4]  that  the  principal 
mechanism  for  the  transfer  of  spectral  energy,  at  least  in  the  limit  of  large  ion 
acoustic  damping  where  free  wave  effects  can  be  ignored,  is  not  resonant  wave-wave 
interactions,  but  rather  the  collapse  of  filaments  which  are  closely  related  to  exact, 
localized,  singular  solutions  of  the  unperturbed  equations.  In  the  two-dimensional 
nonlinear  Schrodinger  equation,  these  filament  solutions  take  the  form  [5] 

where  the  shape”  R  is  the  unique  solution  of  R"  +  tj^R1  —  R  +  R3  =  0,  R'(0)  = 
0,J2(oo)  =  0  without  zeros  in  (0,oo).  The  function  A (t)  approaches  (lnln(i0  - 
—  t)?  as  t  — >  to  [6].  The  filaments  preserve  the  Li  norm  and  so  the  power 
,7*0  =  f  \ip\2dx  carried  by  one  of  these  coherent  structures  is  p  =  2n  R2rjdr] 
and  is  a  pure  number  equal  to  .29  in  our  units  independent  of  time  and  initial 
conditions.  The  value  of  H  for  the  filament  structure  is  exactly  zero  but  each  of  its 
two  components  f  |  V  ^\2dx  and  f  \-ip\*dx  becomes  unbounded.  The  structure  of 
the  singular  filaments  for  the  Zakharov  equations  and  for  the  nonlinear  Schrodinger 
equation  in  three  dimensions  is  somewhat  different  but  in  each  case  the  role  is  the 
same.  These  solutions  are  initiated  by  instabilities  of  the  unperturbed  equations 
and  serve  to  carry  the  energy  from  large  to  small  scales  simply  by  squashing  it  in 
collapsing  filaments.  In  reality,  the  singularity  is  never  reached.  Once  the  filament 
diameter  is  of  the  order  of  the  Debye  radius,  Landau  damping  becomes  important 
and  the  energy  is  transferred  from  the  wave  field  to  electrons,  a  process  known  as 
burnout  This  cycle  of  events  is  shown  in  Figure  1. 

The  instabilitk  •  which  give  rise  to  collapsing  filaments  in  (3)  are  closely  related 
to  the  naturally  occurring  instabilities  of  the  unperturbed  nonlinear  Schrodinger 
equation,  an  infinite  dimensional,  Hamiitonian  system  with  Hamiltonian  H  and  the 
additional  motion  constants  n 0  and  J  —  i  f(rj;  v  ip*  -  ip*  y  ip)dx.  Although  we 
will  discuss  the  geometry  of  phase  space  in  more  detail  in  the  next  section,  it  is 
important  to  mention  at  least  two  types  of  unstable  critical  points  of  the  dynamical 
system.  The  first  are  nonlocal  solutions  which  are  constant  or  purely  periodic  in 
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space  with  a  time  behavior  depending  on  amplitude  (so  they  are  really  unstable  pe¬ 
riodic  orbits).  These  states  axe  unstable  to  sideband  disturbances  (fluctuations  with 
a  spatial  period  close  to  that  of  the  unstable  solutions),  an  instability  known  as  the 
modulaiional  instability.  The  second  type  are  called  cavitons  and  are  analogous  to 
the  solitons  of  the  one  dimensional  equation.  Like  the  collapsing  filament,  they  have 
zero  H  value  and  any  perturbation  which  decreases  H  will  render  them  unstable. 
Because  of  the  existence  of  unstable  solutions,  the  phase  space  of  the  unperturbed 
Hamiltonian  system  contains  many  saddle  points.  Increasing  the  size  of  the  do¬ 
main  or  the  amplitude  of  the  xj)  field  (the  important  parameter  is  the  Li  norm)  will 
create  saddle-center  bifurcations  and  increase  the  number  of  saddles.  Hamiltonian 
perturbations  such  as  the  addition  of  a  term  V(x)\xjx\2  to  H  can  create  additional 
island  chains  of  hyperbolic  saddles  and  elliptic  centers  and  in  particular,  if  V(x)  is 
a  localized  potential  well,  these  new  critical  points  can  correspond  to  the  just  men¬ 
tioned  localized  caviton  states  [7],  The  presence  of  a  small  amount  of  dissipation 
will  turn  centers  into  sinks,  but  the  saddles  will  persist.  Other  non-Hamiltonian  per¬ 
turbations,  such  as  forcing,  can  destroy  the  stable  sinks  and  destabilize  the  weakly 
stable  cavitons  which  axe  localized  in  potential  wells.  In  practice,  this  scenario  of¬ 
ten  arises  [8],  particularly  for  the  Zakharov  equations  and  in  some  instances  for  the 
nonlinear  Schrodinger  equation.  If  ion  damping  is  large,  the  hole  in  the  ion  acoustic 
density  field  remains  after  the  burnout  of  the  electric  field.  This  hole  can  serve  to 
focus  electric  field  energy  and  nucleate  a  metastable  caviton.  However,  as  the  hole 
slowly  collapses,  the  caviton  destabilizes  to  a  collapsing  filament.  We  call  this  the 
nucleational  instability.  It  is  from  a  combination  of  modulational  and  nucleational 
instabilities  that  the  collapsing  filaments  are  bom. 

We  now  turn  to  the  numerical  evidence  for  the  dominance  of  collapse  events. 
Figure  2  displays' the  instantaneous  rate  of  dissipation  7 (t)  =  2^7e(/j)|t/>^.|2,  its 

it 

integral  T(t)  =  —  f  dt'j(t'),  and  the  global  maximum  of  \xf)\2  over  a  time  interval 
including  several  collapses  in  a  strong-turbulent  regime.  —  r  is  just  the  energy  lost 
due  to  dissipation  as  a  function  of  time.  Observe  that  a  large  fraction  of  the  energy 
dissipation  occurs  in  sudden  jumps  and  is  directly  correlated  to  the  collapse  events. 
Energy  is  considered  to  be  dissipated  due  to  collapse  only  if  7  >  70, 

k 

where  70  is  a  threshold  parameter.  If  this  inequality  holds  over  the  interval  [tj,  £2], 
then  the  energy  burnt  out  isT(ti)  —  T(t2),  and  the  total  energy  lost  due  to  collapse 
for  the  entire  simulation  is  the  sum  of  all  such  differences.  Denote  this  total  energy 
lost  due  to  collapse  by  Snc  and  the  total  energy  dissipated  by  Sri?-  In  Figure  3  we 
have  plotted  the  ratio  8nc/6riT  as  a  function  of  70  for  several  different  turbulent 
regimes.  As  70  decreases  the  ratio  increases  and  is  of  course  one  at  70  =  0.  70  =  1 
was  the  smallest  value  that  gave  ratios  judged  to  be  free  from  contamination  by  non¬ 
collapse  events  for  all  cases  studied,  so  it  was  chosen  as  the  standard  cut-off  rate. 
With  this  choice  of  70 ,  as  much  as  80%  of  the  total  energy  dissipated  is  dissipated 
by  collapsing  filaments  in  the  most  energetic  case  examined  (<  n ^  >=  10).  In 
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the  weakest  case  (<  >=  5)  more  than  70%  of  the  dissipated  energy  was  lost  in 

burned  out  filaments.  Furthermore,  these  estimates  of  the  energy  dissipation  are 
conservative  because  we  have  not  yet  included  the  more  gradual  loss  due  to  the  decay 
of  the  high  k  remnant  left  over  from  the  collapse.  The  reason  for  this  remnant  is 
peculiar  to  the  two  dimensional  nonlinear  Schrodinger  equation  and  occurs  because 
a  singular  filament  of  this  equation  carries  the  minimum  threshold  energy  p  required 
to  sustain  collapse.  (In  contrast,  the  energy  of  a  collapsing  Zakharov  filament,  again 
a  constant  of  the  motion,  cam  take  on  a  continuous  range  of  values  above  the  critical 
threshold.)  Therefore,  when  a  significant  portion  of  the  filament’s  spectral  energy 
lies  in  the  dissipative  range  (fc  >  kj),  some  of  its  energy  is  lost,  collapse  is  arrested 
and  only  a  partial  burnout  occurs.  A  histogram  showing  the  distribution  of  the 
fraction  of  the  filament  energy  lost  in  the  burnout  is  shown  in  Figure  4.  It  leaves 
behind  a  remnant  in  the  form  of  broadened  concentric  cylindrical  shells  of  field 
energy  centered  at  the  collapse  site  (see  Figures  le,  f).  In  order  to  follow  the  energy 
associated  with  the  remnant,  we  monitor  (Figure  5)  the  field  energj  inside  a  small 
cell  centered  at  the  collapse  site.  We  note  in  Figure  5b  that  the  burnout  (t  ~  .OS) 
of  the  central  part  of  the  filament,  or  core,  causes  the  fastest  depletion  of  energy 
in  the  cell,  but  also  that  it  is  immediately  followed  by  a  slower  depletion  (the  drop 
over  the  time  interval  .08  <  t  <  .12,  which  we  call  the  “shoulder”).  The  slower  loss 
is  primarily  due  to  the  fluxing  of  the  remnant  through  the  cell  boundary.  If  we  add 
the  energy  loss  rp,  0  <  r  <  1,  due  to  the  burnout  of  the  core  to  that  carried  out  of 
the  cell  by  the  remnant,  we  obtain  a  total  loss  of  energy  in  the  cell  very  close  to  the 
total  amount  p  carried  by  the  collapsing  filament. 

The  remnant  plays  two  roles  in  the  dissipation  process.  First,  it  provides  a  nu¬ 
cleating  center  for  new  collapses,  a  fact  we  have  verified  by  examining  the  spatial 
distribution  of  collapse  sites.  As  a  consequence  the  frequency  of  collapses  increases. 
In  Figure  6,  we  plot  a  histogram  of  the  times  between  successive  eventsw  for  the 
case  76  =  2.  The  distribution  is  nearly  Poisson  with  a  mean  time  between  events 
of  w"1  =<  2  >~  .08.  Second,  the  gradual  damping  of  the  remnant  enhances  the 
ambient  dissipation  rate  <  7,4  >  (i.e.,  the  dissipation  rate  averaged  over  time  inter¬ 
vals  free  from  collapse)  over  that  amount  we  would  expect  from  other  dissipative 
mechanisms.  In  order  to  verify  this,  we  carried  out  the  following  relaxation  ex¬ 
periment.  By  changing  the  sign  of  the  nonlinearity  in  equation  (3)  from  plus  to 
minus,  we  obtain  a  system  with  the  same  linear  dispersion  as  the  plus  system  and 
therefore  presumably  the  same  resonant  wave  interaction  character  but  without 
the  unstable  modulational  and  nucleational  instabilities  and  collapsing  solutions. 
Initializing  the  two  systems  identical!  ,  we  ran  them  with  the  forcing  turned  off 
(76  =  0)  and  compared  their  dissipation  rates  between  collapses  of  the  (+)  system, 
before  their  energies  had  separated  by  more  than  10%  of  the  common  intial  value. 
The  ambient  dissipation  rate  of  the  (+)  system  was  typically  50%  greater  than  that 
of  the  non-collapsing  (  )  system  while  the  collapse  remnants  were  dissipated.  In¬ 
cluding  collapse  events  in  the  tally,  the  (+)  system  lost  energy  three  times  as  fast 
as  the  (-)  system.  We  conclude,  therefore,  that  the  amount  of  energy  dissipated  by 


7 


non-coliapse  mechanisms  is  at  most  one  quarter  of  the  total  energy  dissipated. 

Therefore,  the  energy  dissipation  budget  is  as  follows.  The  average  total  dissipa¬ 
tion  rate  <  7  >  is  given  by 

<  7  >=  ujc  <  r  >  p+  <  >  .  (5) 


The  ambient  dissipation  rate  can  be  eliminated  in  terms  of 
and  we  obtain 

OTl'p 

<  7  >—  u )cp  <r>  ~ — 

OTIq 


<r>P  _ 


In  (5)  and  (6),  <  7  >  depends  only  on  the  mean  turbulence  level  <n^  >  but,  as 
we  have  pointed  out,  <  7.4  >  and  uc  depend  both  on.  <  n,p  >  and  the  nature  of 
the  dissipation  <  r  >  depends  only  on  the  latter.  In  Figure  7,  we  plot  both  <  7  > 
and  u>c  as  a  fraction  of  <  >  and,  consistant  with  (6),  these  curves  are  parallel 

over  a  large  range  of  <  Indeed,  we  checked  the  product  <  7  >  for 

a  range  of  values  of  <  >  and  found  it  to  be  almost  independent  of  the  mean 

turbulence  energy  level  and  equal  to  13  and  equal  to  the  product  <  r  >  taken 

from  Figures  3  and  *  (<  r  >-  'S’  -  .75).  We  expect  that  as  the  dissipation 

mechanism  is  postponed  to  higher  k,  the  fraction  of  energy  lost  in  burnout  will 
increase  and  the  histogram  in  Figure  4  will  accumulate  at  p  =  .29.  We  also  expect 
that  uc  will  decrease  proportionately  and  that  in  this  limit  the  prod*  “'.t  uic  <  r  > 
can  in  principal  be  activated  theoretically  by  arguing  that  in  a  homoclinic  excursion 
the  phase  point  spends  the  longest  time  in  its  cycle  near  the  saddles  in  the  finite 
part  of  phase  space.  If  this  were  the  case  u >c  <  r  >  would  be  given  by  a  product 
of  cr,  the  growth  rate  of  the  modulational  or  nucleational  instabilities,  and  k2,  the 
mean  density  of  collapses.  Both  these  quantities  should  depend  mainly  on  <n^>. 
Verification  of  these  suggestions  requires  extensive  computation  which  is  presently 
underway  but  not  yet  complete. 

In  summary,  then,  we  have  demonstrated  that  the  turbulent  transport  of  field 
energy  to  dissipative  spatial  scales  in  the  two  dimensional  nonlinear  Schrodinger 
equation  is  overwhelmingly  dominated  by  the  explosive  collapse  of  localized  states 
or  filaments,  just  as  it  is  dominated  by  caviton  collapse  in  the  Zakharov  equations. 
To  help  visualize  these  events  as  fluctuations  about  a  basic  homoclinic  excursion(s) 
we  have  plotted  flow  trajectories  projected  onto  the  global  observables  n^,H  and 
l^iir.cs  produce  a  curve  in  three  dimensions  shown  in  Figure  8a.  In  Figure  8b, 
we  have  projected  this  curve  onto  the  three  coordinate  planes  to  aid  visualization. 
Those  trajectories  which  are  fluctuations  about  homoclinic.  excursions  are  the  sparse 
loops  passing  through  the  largest  values  of  |^|£,az.  The  denser,  more  jittery  part  of 
the  curve  lies  in  what  we  call  (in  the  next  section)  the  hash  part  of  the  attractor. 

Whereas  these  results  confirm  the  thesis  that  coherent  collapse  events  dominate 
dissipation,  the  expense  involved  in  two  dimensional  experiments  did  not  allow  us 
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to  run  a  sufficient  number  of  simulations  in  order  to  determine  in  what  limit,  if 
any,  approaches  unity.  In  connection  with  this  question,  and  in  parallel  with 
the  ideas  of  DiPema  and  Majda  [9],  one  would  also  like  to  develop  the  notion  of 
a  weak  solution  in  which  solutions  of  the  undamped  equations  could  be  continued 
beyond  the  collapse  time  by  simply  deleting  the  collapsed  filament  and  lowering  the 
L2  norm  of  the  solution  by  a  fixed  amount.  In  order  to  examine  these  questions, 
we  simulated  the  one  dimensional  nonlinear  Schrodinger  and  Zakharov  models  with 
quintic  nonlinearities  on  a  grid  of  1024  points.  Aliasing  errors  were  removed  by 
smoothly  interpolating  ip  onto  a  grid  of  8  X 1024  points  before  forming  the  nonlinear 
frequency  |i/>|4.  The  results  for  NLS  showed  that,  for  the  same  turbulence  levels, 
the  percentage  of  the  dissipation  rate  accounted  for  by  collapses  rose  from  59%  at 
ha  =  128  (where,  because  the  dissipation  effects  are  clearly  felt  in  the  early  stages, 
there  are  many  failed  attempts  to  form  collapsing  filaments),  to  72%  for  kd  ==  256, 
to  82%  for  kd  =  512.  The  average  loss  of  energy  per  collapse  decreases  with  kd  (the 
amount  of  energy  greater  than  kd  in  the  Fourier  transform  of  the  collapsing  filament 
decreases  with  increasing  kd)  but  the  frequency  of  events  increases  proportionately 
so  that  the  average  dissipation  rate  remains  the  same.  The  evidence  clearly  suggests 
that  in  the  large  kd  limit,  sill  the  energy  dissipated  is  dissipated  by  collapse  events 
which  occur  infinitely  often  with  infinitisemal  losses  of  energy  per  event. 

A  much  more  striking  result  was  obtained  when  we  ran  the  very  same  experiment 
on  the  Zakharov  model 

ipt  -  i'Pxx  +  ip*P  =  F-eD  (7) 

p„  +  2z/opt-pIX  =  (|V>|4)XI  (8) 


where  F  and  D  are  as  before  and  2v  o  pt  is  a  convolution  integral  modelling  ion 
damping.  The  ion  acoustic  field  p(x,f)  is  no  longer  slaved  to  the  electric  field 
intensity  (for  NLS,  p  in  the  equation  (7)  is  replaced  by  — |r/»|4).  Therefore,  during 
collapse,  in  which  the  fields  take  on  a  self-similar  form  close  to  (2)  but  in  which 
the  inertial  acceleration  p»  is  also  important,  the  cavity  formed  by  the  ion  field 
encourages  total  burnout  of  the  filament.  We  found  that  when  e  =  1,  kd  =  32,  < 
ny  >=  1.625,  the  average  energy  <r  >  p  lost  per  event  was  1.95.  the  average  time 
between  events  was  .83  and  ^  was  92%.  (For  the  Zakharov  model,  the  amount 
of  energy  carried  in  the  filament,  a  constant  of  the  motion  for  the  unperturbed 
equations,  can  take  on  a  range  of  values  greater  than  the  threshold  value  p  =  .43  in 
this  case.  Thus  r  >  1).  The  distribution  of  r  rose  sharply  after  r  =  1,  had  a  mean 
of  1.95/.43  and  had  a  relatively  long  tail.  Also,  we  observed  that  all  the  collapses 


red  in  9  nonoverlapping  sites,  were  driven  by  the  nucleational  instability  and 


the  sites  drifted  about  the  box.  When  we  increased  kd  to  64,  <  >=  1.8,  the 

frequency  of  events  increased  to  (.66)_1  (there  were  151  events  in  100  time  units) 


and  the  distribution  of  r  came  closer  to  one.  The  average  energy  lost  per  event  was 
1.55  and  was  95%.  (In  contrast,  at  the  same  parameter  values  and  in  the  same 
time  interval,  the  NLS  model  had  972  events  with  an  average  energy  loss  of  .09  and 
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was  72%. )  Further,  when  we  decreased  e  to  j  (again  kj  =  64),  the  distribution 

of  r  came  closer  to  threshold  and  jjf-  was  9S%.  This  leads  us  to  conjecture  that 
as  kd  increases  and  c  decreases,  the  distribution  of  r  will  cluster  about  one,  and 
that  each  collapsing  event  will  bum  off  exactly  the  threshold  energy  p  =  .43.  In 
this  asymptotic  limit,  the  weak  solution  of  the  unforced,  undamped  equation  in 
which  the  turbulence  eventually  decays,  is  found  by  simply  removing  the  collapsing 
filament  from  the  field  and  reducing  the  L2  norm  by  p.  (The  weak  limit,  i.e.  ?/>(x,  0) 
is  a  weak  solution  if  J  ip(x)ip(x1 0 )dx  —  lim  f  (p(x)Tp(x ,  A )dx  where  (p  is  smooth,  of 

a  collapsing  filament  is  zero  because  it  oscillates  very  fast  and  its  width  decays  at 
a  faster  rate  than  its  amplitude  increases).  We  will  report  more  details  elsewhere. 
We  are  also  currently  testing  some  ideas  concerning  the  estimation  of  the  frequency 
of  events. 

In  Figure  9,  we  give  the  angle- averaged  correlation  function 

f 2ir  j/\ 

Gi>(p)  =  Jo  2^r  <  (*  +  W(*)  > 

= jf*  <  hW)i2  > 

k 

=  £ MkP)  <  |^(S)p  > 

k 

in  the  turbulent  regime  75  =  .2.  Observe  that  the  correlation  length  A  =  2 ir  < 
k2  .7  is  much  less  than  the  simulation  box  size  (2tt)  so  that  the  turbulence 

is  independent  of  the  periodic  boundary  conditions. 

We  end  this  section  with  a  few  remarks  about  the  two-dimensional  numerical 
scheme.  In  calculating  solutions  to  (3),  we  use  the  following  “split-step”  numerical 
algorithm  to  advance  the  solution  from  t  to  t  +  dt. 

i()(x,t  +  dt)  =  ^^1^(x)exp(-i|^>^1^(x)j2di) 

where  (*  refers  to  Fourier  transform) 

=  1 p(k,  t)exp(~fb(k)  -  7e(fc)  +  ik2)dt 


This  method  exactly  conserves  n^(=  f  dx\ip\2)  in  the  conservative  limit  (7e  -+  0) 
and  when  aliasing  errors  are  absent;  errors  in  the  Hamiltonian  =  f  (j  V0|2 -y|^|4)dx 
are  0 (dt2).  Aliasing  errors  incurred  during  the  nonlinear  part  of  the  evolution  are 
minimized  by  strong  dissipation  at  short  scales  (see  figure  9)  that  curtails  the  time- 
averaged  spectrum  (see  figure  10).  Our  time  step  (dt)  is  10“4,  and  is  sufficient  to 
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resolve  the  maximum  nonlinear  frequencies  {\ip\2max  <  4  x  103)  encountered  in  the 
most  intense  turbulent  regime  studied.  The  Landau  damping  function  used  is  shown 
in  Figure  10.  Its  asymptotic  shape  is  chosen  to  be  sufficiently  steep  (~  &1+£,  e  >  0) 
so  that  the  collapsing  filament  is  eventually  arrested  but  not  too  steep  so  as  to  reflect 
the  collapsing  filament  in  k  —  space.  It  depends  on  a  single  parameter,  the  ion-to- 
electron  mass  ratio  M,  which  in  our  analysis  we  take  to  be  7344  (the  mass  ratio  in 
a  singly-ionized  Helium  plasma),  and  increases  abruptly  at  k  =  .2 k&  —  .2(| \/M). 
The  turbulence  energy  level  is  changed  by  varying  the  forcing  term  F.  In  all  cases 

F  =  jb[ifi(Akx,  Aky )  +  ip(-  A  kXl  A ky)  +  tp{A kx,  -  Aky)  +  rp(-  A  kX)  -  A  fcj,)], 


i.e.,  energy  is  injected  by  linearly  destabilizing  four  modes  symmetrically  in  k  — 
space ,  where  A kx  =  A ky  =  and  Lx  =  Ly  =  2ir.  (Our  simulations  are  in  a 
square  box  2ir  long  on  each  side  spanned  by  a  regular  grid  of  1282  points.)  Because 
this  scalar  forcing  is  isotropic,  the  turbulence  is  isotropic.  (See  Figure  10).  Notice 
that  as  a  filament  collapses  it  will  decouple  from  the  forcing  which  is  confined  to  long 
wavelengths.  Thus  the  dominant  transport  mechanism  (collapse)  is  independent  of 
the  details  of  injection. 

Geometry  of  phase  space 


While  the  full  attractor  A  is  large  dimensional,  the  dominance  of  filaments  sug¬ 
gests  that  A  is  the  fuzzy  covering  of  a  skeleton  which  consists  of  generalized  saddle 
points  M  and  S  and  two  heteroclinic  orbits.  In  the  case  study  just  discussed,  M 
is  the  set  of  saddle  points  which  are  modulationally  or  nucleationally  unstable  and 
S  is  the  idealized  collapsing  filament,  namely  the  singular  self-similar  solution  of 
the  two  dimensional  nonlinear  Schrodinger  equation.  The  unstable  manifold  of  M 
intersects  the  stable  manifold  of  S  and  vice  versa.  Based  on  this  picture,  A  can  be 
subdivided  into  two  subsets  Ah  and  Ac .  When  the  system  is  in  the  turbulent  soup 
Ah,  it  is  dominated  by  what  we  call  the  hash  modes  consisting  of  the  background 
field  and  radiation  modes  left  over  from  the  formation  of  the  coherent  structures 
(the  filaments  of  Langmuir  turbulence).  This  set  is  generally  large  dimensional  and 
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in  the  case  of  shear  flow  turbulence  could  have  dimension  R t.  On  the  other  hand, 
when  the  system  is  in  Ac,  it  is  dominated  by  the  coherent  structures, . although 
large  dimensional  fluctuations  about  these  orbits  are  still  present.  The  orbits  in 
Ac  are  organized  by  a  low  dimensional  submanifold  of  solutions  B,  which  is  the 
stable  manifold  of  an  unstable  saddle  point  S,  possibly  at  infinity,  corresponding 
to  an  idealized  state.  In  general,  S  may  be  a  manifold  of  idealized  states  and  in 
many  cases  will  be  the  orbit  of  some  group  action  representing  the  symmetries 
(translation,  scaling)  which  are  properties  of  the  system  at  hand.  In  particular, 
if  the  spatial  dimension  is  large,  coherent  structures,  each  with  the  same  shape, 
can  be  localized  at  several  locations,  although  occurring  at  different  times,  with 
an  approximately  uniform  density  reflecting  the  amount  of  power  (energy  density 
times  area  or  volume)  needed  to  form  the  structure.  (In  Langmuir  turbulence,  each 
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filament  requires  an  area  of  k~2\  k  depends  on  the  turbulence  intensity  level.)  In 
constructing  the  phase  space,  one  should  think  of  the  system  as  being  confined  to  a 
box  of  a  sufficient  size  to  contain  one  coherent  structure  or  to  a  lattice  in  which  case 
the  idealized  solution  S  consists  of  an  array  of  discrete  translates  of  the  prototype 
and  we  consider  the  system  to  be  in  the  state  S  when  any  one  is  triggered.  The 
unstable  manifold  of  S  (the  burnout  due  to  Landau  damping)  is  asymptotic  to  Ah- 
The  speeds  of  attraction  and  repulsion  at  S  are  governed  by  two  entirely  different 
processes  and  are  not  related.  The  former  is  determined  by  the  faster  than  expo¬ 
nential  rate  at  which  the  filament  becomes  singular,  the  latter  is  governed  governed 
by  dissipation.  The  picture  we  have  described  is  drawn  schematically  in  Figure  12, 
and  corresponds  closely  to  an  actual  orbit  of  the  nonlinear  Schrodinger  equation 
shown  already  in  Figure  8. 

We  now  turn  to  a  discussion  of  the  structure  of  the  phase  space,  the  nature  of  the 
saddle  point  M  and  its  unstable  manifold.  Whereas  the  perturbations  of  forcing 
and  damping  modify  the  phase  space  structure,  it  is  the  topology  of  the  phase  space 
of  the  unperturbed  system  which  sets  the  stage  for  the  large  homoclinic  excursions. 
Two  properties  of  the  unperturbed  system,  which  is  Hamiltonian,  are  crucial.  One  is 
the  fact  that  its  phase  space  contains  many  saddle  points  corresponding  to  unstable 
fixed  points  and  periodic  orbits  and  separatrices  joining  these  saddles.  The  second 
crucial  property  is  the  non-compactness  of  the  constant  energy  surfaces  and  the 
existence  of  heteroclinic  connections  (HCI’s)  which  join  certain  of  the  saddle  points 
to  infinity.  Whereas  Hamiltonian  perturbations  of  this  system  can  destroy  resonant 
KAM  tori  and  create  new  chains  of  elliptic  centers  and  hyperbolic  saddles  and  non- 
Hamiltonian  perturbations  can  turn  the  centers  into  attracting  sinks  or  destroy  them 
altogether,  the  hyperbolic  nature  of  the  original  saddles  and  their  HCI’s  remain 
intact  and  play  a  dominant  role  in  the  dynamics  of  the  perturbed  system.  As  an 
example,  consider  the  perturbed  Duffing’s  oscillator 

x"  +  x  —  2x3  =  F(t)  —  (—  )nx'  (9) 

Xo 

where  F(t)  is  small  and  the  additional  factor  (~)n  has  been  added  to  the  damping 
term  in  order  to  arrange  for  it  to  turn  on  only  when  x  is  large.  The  unperturbed, 
exactly  integrable  Hamiltonian  system  represents  a  particle  in  a  quartic  potential 
V(x)  =  |(x2  —  x4)  with  critical  points  x1  =  x  =  0  (a  center)  and  x‘  —  0,  x  =1 
(saddle  points).  Its  phase  plane  is  shown  in  Figure  13.  The  unstable  and  stable 
manifolds  emanating  from  x  =1  are  the  non-compact  energy  level  surfaces 

E  =  --  =  \x'2  +  |x2  —  \x*  which  join  the  saddles  to  infinity  through  the  HCI’s 
x'  =  (x2  —  |)  for  x  >  0,x  >  and  x  =  — (x2  —  \)  for  x  <  0,x  <  —  The  other 

branches,  the  stable  manifolds  of  x-  and  x  =1  are  repellers  for  positive  time. 

In  particular,  observe  that  all  energy  level  surfaces  E  >  \  and  those  for  E  <  \  for 
which  x2  >  ^  approach  the  HCI’s  at  infinity  and  that  any  point  on  these  trajectories 
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will  reach,  infinity  in  a  finite  time.  This  can  be  seen  in  severed  ways.  First,  notice 
that  the  distance  x\  —  x 2  between  two  constant  energy  curves  E\  and  E2  at  a  given 
value  of  x'  is  given  by  {x\  -  x^)(xf  +  x\  -  1)  =  2E2  -  2E\  so  that  asymptotically 
Xi  —  x2  ~  (2E2  —2Ex)(x\  +X2)-1(xf  +X2)-1.  Second,  one  can  construct  a  Laurent 
series  for  any  solution  in  the  neighborhood  of  infinity, 


which,  because  the  system  is  integrable,  has  the  Painleve  property  (the  series  is 
Laurent  and  has  the  required  number  of  free  constants,  namely  two,  to  and  c; 
all  the  other  coefficients  an  can  be  calculated  explicitly  in  terms  of  c).  The  free 
constant  c  is  related  to  the  energy  E  =  — oc  +  The  distance  between  any  two 
orbits  is  the  minimum  over  t\  and  t2  of  x(fi,ci)  -  x(t2,c2 )  which  goes  to  zero 
like  (f  —  fi)3.  We  have  thus  introduced  a  new  concept  into  Hamiltonian  systems 
with  non-compact  energy  surfaces ,  namely  the  notion  that  certain  orbits  (HCI’s) 
can  serve  as  attractors  for  all  other  orbits  at  infinity  in  the  sense  that  they  are 
asympt  .ti'cally  all  the  same.  This  notion  does  not  violate  conservation  of  volume  in 
phase  space.  Close  to  infinity,  the  stretching  of  nearby  points  in  a  direction  parallel 
to  the  HCI  is  faster  than  exponential  (the  local  stretching  exponent  in  this  direction 
is  infinite)  and  therefore  the  contraction  in  the  direction  perpendicular  to  the  HCI  is 
also  faster  than  exponential.  For  small  forcing,  the  orbits  of  the  perturbed  systems 
which  escape  to  infinity  also  follow  the  HCI  very  closely  until  x  >  x0  at  which  point 
dissipation  sets  in.  However,  the  important  point  is  that  for  x0  large,  all  escaping 
orbits  are  approximately  parallel  to  the  HCI  over  large  regions  of  the  phase  space. 

One  can  look  at  the  nature  of  the  phase  plane  at  the  line  at  infinity  by  introducing 
homogeneous  projective  coordinates  given  by  x  =  — ,x'  =  y  =  j  in  which  coordi¬ 
nates  the  direction  field  (x-2x3)dx+yciy  =  0  becomes  (X  Z3 —2X*  Z)dX +Y  Z3  dY + 
(2X4  —  X2Z2  —  Y2Z2)dZ  =  0.  Since  the  critical  point  at  infinity  is  X  =  Z  =  0, 
it  is  interesting  to  look  at  the  direction  field  in  the  affine  chart  given  by  Y  ^  0. 
The  line  field  in  the  X ,  Z  coordinates  is  obtained  by  setting  Y  —  1  and  dY  =  0. 
We  observe  in  Figure  12  that  invariant  manifolds  ( Y 2  +  X2)Z2  —  X 4  —  2EZ4  =  0 
organize  the  flow  field  in  the  neighborhood  of  the  line  at  infinity  and  in  particular 
all  of  them  converge  from  the  half  plane  X  >  0  to  the  point  X  —  Z  =  0  on  the 
line  at  infinity  and  then  reemerge  in  the  left  hand  plane.  The  convergence  of  these 
curves  at  infinity  again  shows  that  the  HCI  is  an  attractor  there. 

This  simple  example  illustrates  how  the  topology  of  the  phase  space  of  the  un¬ 
perturbed  system  is  important  in  controlling  ,.he  dynamics  of  the  perturbed  sys¬ 
tem.  Likewise  in  the  weakly  perturbed  nonlinear  Schrbdinger  equation,  in  which 
forcing  is  applied  at  small  wavenumbers  and  the  damping  at  large,  it  is  the  topol¬ 
ogy  of  the  phase  space  struct  ure  of  the  unperturbed  problem  which  dominates  the 
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dynamics  of  the  perturbed  problem.  The  unperturbed  problem  is  again  Hamil¬ 
tonian  although,  in  this  case,  it  is  infinite  dimensional  with  motion  constraints 
rty  =  / ipiJ)*dx,J  =  i  ftp's?  rb'-^'Y^dx  and  H  -  \  /(vV’1  VV  -  %tp2tp*2)dx.  In 
particular,  as  we  have  mentioned,  it  is  known  that  the  energy  surface  is  non-compact 
and  when  H  <  0,  there  are  orbits  in  any  surface  n^,  >  nc  =  2tc  J0°°  R2(rj)T]drj  which 
connect  the  unstable  saddle  points  to  infinity.  We  remark,  however,  that  in  our 
simlulation  H  was  almost  always  positive  although  taken  in  the  neighborhood  of  a 
collapse  it  was  zero.  The  negativity  of  H  is  a  sufficient  but  not  necessary  condition 
for  collapse. 

Our  picture  of  the  phase  space  then  is  as  follows.  For  low  levels  of  forcing  (either 
F  is  small  or  the  domain  area  is  small),  the  modulus  of  the  electric  field  grows 
till  >  nc  at  which  stage  the  phase  point  eventually  comes  close  to  a  modula- 
tionally  unstable  saddle  point  which  is  joined  to  infinity  via  a  collapsing  filament. 
The  instability  sets  in,  the  filament  is  formed,  it  collapses,  a  certain  portion  of  its 
energy  is  dissipated  (depending  on  the  damping  structure)  and  a  remnant  of  high 
k  waves  is  left  over.  The  system  returns  close  to  its  original  state  and  the  process 
is  approximately  repeated  with  the  initialization  of  the  collapsing  filament  being 
due  to  either  modulational  or  nucleational  instabilities.  In  particular,  the  remnant 
can  provide  a  cavity  for  the  nucleation  of  metastable  cavitons  which  collapse  once 
the  support  disperses  away.  For  larger  values  of  the  applied  stress  (either  larger 
forcing  or  larger  boxes),  the  value  n  of  the  mean  turbulence  level  is  much  larger 
and  can  be  many  times  nc.  In  this  case,  spatial  correlations  decay  rapidly  and 
many  collapsing  filaments  can  occur  at  different  spatial  locations  although  since 
each  event  is  so  rapid,  they  will  rarely  occur  in  the  same  time  intervals.  For  these 
large  turbulence  levels,  the  saddle  point  M  in  Figure  12  represents  a  collection  of 
saddle  points  in  the  turbulent  soup  part  of  the  attractor  Ah,  some  of  which  cor¬ 
respond  to  analogues  of  the  modulational  instability  of  more  complicated  periodic 
shapes  and  others  of  which  are  best  understood  as  nucleational  instabilities.  As 
we  have  mentioned,  the  collapsing  filament  is  often  initiated  in  the  remnant  of  a 
hole  left  by  a  previous  event  or  in  holes  simply  produced  by  random  fluctuations. 
No  matter  what  the  mechanism,  the  important  point  is  that  the  attractor  contains 
many  saddle  points  M  which  are  connected  to  infinity  through  HCI’s  and  the  phase 
point  passes  sufficiently  close  to  these  points  so  that  the  large  homoclinic  excursion 
is  initiated. 

We  mention  here  that  a  similar  situation  obtains  in  the  Euler  equations.  For 
low  values  of  the  applied  stress,  laminar  states  can  be  destabilized  by  identifiable 
instability  mechanisms,  centrifugal  instabilities,  mean  flow  profiles  with  inflexional 
points  and  so  on.  These  instabilities  do  not  go  away  when  the  fluid  becomes  fully 
turbulent.  Indeed  saddle  points  representing  local  inflexional  and  centrifugal  in¬ 
stabilities  remain  very  much  part  of  the  strange  attractor  of  the  high  Reynolds 
number  Navier  Stokes  equations  and  play  a  large  role  in  transferring  energy  to  high 
wavenumbers  where  dissipation  acts.  In  highly  turbulent  flows,  the  inflexional  pro- 
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file  of  the  mean  flow  is  not  sustained  for  all  time  uniformly  in  space,  but  if  it  is 
sustained  long  enough  in  a  local  region  (see  the  discussion  on  turbulent  boundary 
layers  in  the  following  section),  then  rapidly  growing  packets  of  three  dimensional 
inflexional  instabilities  can  erupt  and  carry  energy  off  to  the  dissipation  cemetaries. 

Other  illustrations  of  the  random  occurrence  of  coherent  events 

This  view  of  turbulent  transport,  in  which  the  major  contribution  comes  from 
a  part  of  the  attractor  which  corresponds  to  coherent  structures  whose  dynamics 
can  be  captured  with  a  finite  number  of  degrees  of  freedom,  while  not  universal, 
should  have  widespread  application.  We  now  suggest  several  other  contexts  in  in 
which  these  ideas  may  be  important.  We  dress  that  in  no  case  do  we  feel  we  have 
presented  the  complete  and  final  solution,  but  simply  wish  to  suggest  that  it  may  be 
worthwhile  to  look  at  each  of  these  situations  from  a  new  viewpoint.  The  first  is  the 
burst-sweep  cycle  [10]  which  appears  to  dominate  the  production  of  turbulent  energy 
and  momentum  transport  in  turbulent  boundary  layers.  Although  no  candidate  for 
the  coherent  structure  which  models  the  four  step  process  constituting  the  burst 
phase  of  this  cycle  has  yet  been  put  forward  (there  have  been  suggestions  [11]),  it  is 
not  improbable  that  a  singular,  finite  time  collapse  solution  of  the  Euler  equations 
plays  a  significant  role.  This  basic  structure  could  probably  be  captured  by  a  simpler 
set  of  p.d.e.’s,  analogous  to  the  Zakharov  equations,  in  which  a  long  scale,  three- 
dimensional  Tollmein-Schlichting  like  disturbance  with  large  downstream  vorticity 
(streaks)  is  synchronously  coupled  to  a  packet  of  short  inflexional  waves  produced 
by  an  inflexional  instability  of  the  mean  turbulence  velocity  profile  caused  by  the 
induced  outflow  of  the  streaks  [10,11].  The  inflexional  instability,  ever  present  in 
shear  flows,  together  with  the  finite  amplitude  distortion  of  the  mean  profile  by  the 
long  scale  state,  plays  the  role  that  modulational  and  nucleational  like  instabilities 
do  in  Langmuir  turbulence.  It  continuously  drives  the  initial  stage  of  the  bursting 
solution. 

Convection  at  large  Rayleigh  numbers  is  a  second  area  in  which  the  transport 
properties  axe  likely  to  be  organized  by  coherent  structures,  in  this  case  the  slightly 
tilted  plumes  observed  by  Krishnamurti  and  Howard  [12]  which  dominate  the  flow 
for  a  large  range  of  Prandtl  P  and  Rayleigh  numbers  Ra-  More  recent  observations 
by  Libchaber  [13]  iD  what  he  calls  the  soft  turbulence  regime  (Rz  <  !07)  would 
seem  to  lend  some  credence  to  che  admittedly  simplified  picture,  suggested  many 

years  ago  by  Howard  [14],  that  the  plumes  would  occur  at  a  frequency  determined 
—  2  _1 
by  the  time  Ra  J  it  takes  for  a  conductive  layer  to  build  to  a  depth  Ra  ?  so  as  to 

be  convectively  unstable.  A  crude  picture  in  which  cne  imagines  that  the  coherent 

plumes  carry  oft  all  the  heat  in  the  unstable  conductive  layer  would  suggest  a  Nusselt 

number  proportional  to  R$  which  is  approximately  what  is  observed.  The  inferences 

of  the  Howard  argument  are  also  consistent  with  a  “mixing-length”  analysis  given  in 

an  earlier  paper  by  Kraichnan  [15].  To  be  sure,  this  picture  is  greatly  oversimplified 

and  it  is  hard  to  see  how  the  plumes,  which  combine  into  larger  and  larger  ones, 
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remember  the  Howard  time  scale.  Nevertheless,  Libchaber  observes  this  scale  as 
dominant  in  the  power  spectrum  of  his  signal,  measured  many  conductive  layer 
depths  from  the  bottom  of  the  tank,  over  a  large  range  of  Rayleigh  numbers  up  to 
approximately  107.  At  this  stage,  krai  Reynolds  numbers,  which  are  proportional 
to  the  square  root  of  the  Rayleig1  _mber,  are  in  the  shear  flow  instability  range 
and  other  structures  are  produce  vhich  modify  somewhat  the  Howard  picture. 

Third  [15],  we  suggest  that  the  self-similar  features  seen  in  the  Rayleigh-Taylor 
unstable  interface  two  fluids  of  different  densities  are  exact,  singular  solutions  of  the 
Euler  equations,  which  if  they  were  known,  would  help  one  to  calculate  an  approxi¬ 
mation  to  the  mass  flux  of. heavier  through  lighter  fluid.  In  the  absence  of  a  length 
scale  (surface  tension  is  ignored),  the  set  of  attracting  states  S  at  infinity  will  be  a 
complicated  fractal  set  reflecting  the  scale  symmetries,  but  the  basic  elements  of  S 
would  represent  the  spikes  (containing  the  heavier  fluid  falling  through  the  lighter) 
and  bubbles  (the  lighter  fluid  pushing  up  through  the  heavier).  Fourth  [16],  it  is  not 
unreasonable  to  argue  that  the  dissipation  rate  of  shear  flow  turbulence  arises  not 
from  the  wave- wave  interaction  familiar  from  Fourier  space  cumulant  descriptions, 
but  instead  to  suggest  it  is  dominated  by  the  formation  and  destruction  of  thin 
vortex  sheets  (the  state  S)  in  which  almost  all  of  the  vorticity  is  concentrated.  We 
conjecture  that  the  initial  formation  of  surfaces  of  vorticity  concentration  follows 
singular  solutions  of  the  Euler  equations  driven  principally  by  inflexional  instabili¬ 
ties.  Once  formed  (the  system  is  close  to  5),  these  sheets  are  notoriously  unstable, 
to  Kelvin-Helmholz  instabilities  of  the  tangential  velocity  discontinuities  and  to 
Taylor- Gortler  centrifugal  instabilities  of  the  helical  flows  induced  by  vortex  lines 
embedded  in  the  curved  sheets.  These  secondary  instabilities  [17,18]  quickly  trans¬ 
fer  the  energy  to  the  viscous  cemetaries  where  energy  is  dissipated- and  the  system 
is  returned  close  to  a  state  where  approximately  the  same  cycle  can  repeat. 

Our  last  case  study  concerns  the  generation  of  turbulent  spots  in  a  boundary 
layer.  This  example  is  somewhat  different  in  nature  to  the  preceding  examples 
but  nevertheless  exemplifies  the  general  idea  of  transport  occurring  as  a  result  of 
the  random  occurrence  of  coherent  events  -  in  this  case  the  formation  and  convec¬ 
tion  of  the  spots  -  which  can  be  associated  with  homoclinic  excursions.  We  shall 
also  discuss  a  simplified  model  and  make  some  experimental  predictions.  However, 
before  doing  this,  let  us  recall  the  initial  stages  of  the  boundary  layer  instability 
without  receptivity  (undue  influence  of  fluctuations  in  the  outer  flow).  When  the 
Reynold’s  number,  based  on  the  boundary  layer  thickness  which  depends  on  the 
distance  from  the  leading  edge,  reaches  a  critical  value,  the  unstable  two  dimen¬ 
sional  Tollmein-Schlichting  (T.S.)  wave  concentrated  at  the  critical  layer  quickly 
becomes  three  dimensional  thereby  creating  downstream  vorticity.  The  induced 
secondary  flow  causes  an  outward  flow  near  the  crests  (the  most  downstream  part 
of  the  distorted  “  vortex”  lines)  and  an  inward  flow  near  the  troughs.  The  three 
dimensional  distortion  is  enhanced  by  the  mean  shear,  the  downstream  vorticity 
increases  and  the  initially  sinusoidal  three  dimensional  deformation  of  the  vortex 
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line  becomes  strongly  distorted  into  the  form  of  an  elongated  hairpin  vortex.  The 
induced  outward  flow  carries  slow  moving  fluid  to  the  outer  part  of  the  boundary 
layer  and  causes  the  mean  turbulent  profile  to  develop  inflexion  points  localized 
in  both  time  and  the  spanwise  direction.  At  these  points,  short  wavelength  pack¬ 
ets  of  an  inflexionally  unstable  nature  can  grow  and  be  absolutely  (rather  than 
convectively)  unstable  if  they  can  phase  lock  with  the  travelling  three  dimensional 
disturbance  which  gives  rise  to  their  formation  in  the  first  place.  Along  these  tra¬ 
jectories  turbulent  spots  are  formed  which  spread  spatially  into  surrounding  regions 
by  destabilizing  the  neighboring  laminar  flow.  Our  picture  of  the  spots  is  that  they 
are  essentially  localized  wavepackets  in  which  the  waves  are  breaking.  Suppose  one 
monitors  the  velocity  field  v(t )  at  a  point  just  before  the  area  in  which  the  spots 
tend  to  form.  A  reconstruction  of  the  attractor  from  this  time-series  should  show 
an  attractor  A  with  a  well-defined  phase  corresponding  to  the  phase  of  the  incom¬ 
ing  T.S.  waves  and  we  propose  using  the  methods  of  the  following  section  to  check 
thiis  conjecture.  However,  monitoring  the  velocity  field  in  the  full  region  one  will 
find  an  attractor  A  which  will  have  an  Ah  somewhat  similar  to  A ,  but  also  with 
a  part  Ac  consisting  of  the  homoclinic  excursions  which  occur  when  a  spot  forms 
and  is  convected  downstream.  When  the  spot  reaches  the  region  where  they  have 
coalesced  or  left  the  channel,  the  system  has  returned  to  Ah- 

The  phase  space  associated  with  this  dynamics  is  somewhat  different  from  that 
proposed  previously  for  collapsing  filaments.  Nevertheless,  while  the  homoclinic 
excursion  will  not  go  to  infinity,  it  will  have  several  features  in  common  with  the 
previous  picture.  We  consider  here  a  very  simple  mathematical  model  to  illustrate 
what  we  have  in  mind.  Consider  a  locally  unstable  attracting  periodic  orbit  7  in 
phase  space.  This  is  an  orbit  which  is  attracting  in  the  usual  sense  that  all  its 
Floquet  multipliers  have  modules  less  than  one,  but  which  at  some  points  x  in  7 
is  actually  repulsive.  An  example  of  such  as  orbit  would  be  given  by  the  periodic 
solution  x  =  0  for  the  time-dependent  system. 

dcr  d  s 

X(x,  s) :  =  a(s)x  +  0(x2),  —  =  1,  (s,  s)  6  R  x  T 1  (10) 

at  at 

where  a(s)  is  a  smooth  function  which  is  —  1  fjr  0  <  s  <  1/2, 1  for  5/8  <  s  <  3/4  and 
nearly  linear  in  between.  Another  more  natural  example  is  given  by  the  following 
flow  in  the  plane  (See  Figure  14).  If  the  periodic  orbit  7  is  close  to  the  saddle  M 
then  7  will  be  unstable  near  M  even  though  7  is  attracting  in  the  usual  sense.  Now 
returning  to  the  general  picture,  suppose  that  such  a  system  is  subject  to  some 
periodic  or  almost-periodic  forcing  or  perturbation.  If  the  perturbation  is  large 
when  the  state  is  near  the  unstable  part  of  7  then  one  will  see  a  relatively  large 
excursion  of  7,  otherwise  the  perturbation  will  be  quickly  damped  out.  This  picture 
is  a  useful  oversimplification  of  the  dynamical  phase-space  structure  associated  with 
spots.  The  attractor  Ah  like  7  is  unstable  in  some  parts  so  that  the  incoming  waves 
perturb  the  state  sufficiently  to  send  it  on  an  excursion  in  Ac-  Of  course,  we  are  not 
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suggesting  that  Ac  and  Ah  are  uncoupled  as  in  the  simple  example:  indeed,  they 
are  part  of  a  single  strange  attractor.  But  these  can  be  thought  of  as  coming  from  a 
common  picture  as  we  now  indicate.  Consider  again  the  situation  shown  in  Figure  . 
If  this  system  is  periodically  forced  then  the  total  system  is  described  by  a  Poincare 
map  corresponding  to  advancing  one  period  of  the  forcing  and  the  above  picture 
can  be  regarded  as  being  associated  with  the  vector  field  X(x,s)  given  in  Equation 
(10)  for  s  frozen.  Suppose  that  the  way  in  which  this  frozen  vector  field  changes 
during  one  cycle  of  the  forcing  is  such  that  a  limit  cycle  gets  pushed  through  the 
stable  manifold  of  the  saddle  M.  Then  homoclinic  orbits  are  formed  in  the  complete 
system  and  this  will  have  am  attractor  consisting  of  a  set  Ah  of  orbits  close  to  the 
cycle  7  and  horseshoes  near  M,  and  homoclinic  excursions  going  from  a  region  close 
to  M  back  to  Ah-  If  this  is  a  reasonable  approximation  of  the  phase  space  then 
we  can  make  a  prediction.  There  should  be  a  correlation  between  the  creation  of 
spots  and  the  phase  on  Ah-  Experimentally  this  can  be  checked  by  reconstructing 
A  using  the  usual  time-series  methods  and  recording  when  spots  occur  in  the  region 
just  downstream  from  the  monitoring  point. 

Identification  of  coherent  structures 

As  the  final  part  of  this  paper,  we  want  to  propose  a  method  of  time  series  analy¬ 
sis  which  can  identify  those  modes  which  contribute  most  to  transport.  These  ideas 
are  based  on  and  are  a  modification  of  the  singular  value  decomposition  (SVD) 
methods  of  Broomhead  and  King  [19]  to  improve  phase-space  constriction  tech¬ 
niques  and  Lumley’s  method  for  identifying  coherent  structures  in  turbulent  flows. 
We  will  elaborate  on  these  ideas  in  more  depth  in  a  future  publication  with  David 
Broomhead  but  here  we  will  sketch  the  main  points.  The  principal  goal  of  the 
method  is  to  extract  organized  structures  which  may  come  and  go  from  what  ap¬ 
pears  to  be  a  sea  of  disorganized  and  statistical  fluctuations. 

Firstly,  let  us  recall  the  basic  ideas  from  the  perspective  of  dynamical  systems. 
We  suppose  that  we  have  a  dynamical  system  acting  as  a  Hilbert  space  if  with 
an  attractor  A  which  has  a  natural  measure  v.  Let  the  mean  v  =  f  vv(dv)  be 
subtracted  from  each  of  the  vectors  v  so  that  they  then  have  zero  mean.  We  define 
an  orthogonal  basis  e\,  e2>  •  •  •  of  if  by  maximizing 

F  =  J  (e  -  v)2v(dv)  +  A(|e|2  —  1)  (11) 

to  get  ej,  and  then,  given  ei,...,e,-,  taking  for  e,+i  the  maximum  of  F  in  the 
orthogonal  complement  of  , . . . ,  e j .  We  ignore  for  the  moment  the  non-generic 
situation  where  F  does  not  have  such  maxima  and  where  A  is  contained  in  a  finite- 
dimensional  linear  subspace.  For  the  maximal  basis,  we  find 

F  =  of  -  J(e{  •  vfv{dv)  (12) 
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Then  the  <7*  can  be  regarded  as  the  moments  of  inertia  of  A  if  A  is  regarded  as  a 
body  the  mass  of  whose  points  is  given  by  v.  The  e,  are  the  corresponding  axes 
of  inertia  and  these  can  be  easily  calculated  by  constructing  the  inertia  tensor  by 
differentiating  (11)  and  solving  for  its  eigenvectors.  In  this  case  the  inertia  tension 
is  the  nxn  matrix  X*X  where  X  is  the  nxN  matrix  whose  ith  row  is  the  vector  vl. 
This  is  completely  practical  procedure  which  we  have  carried  out  on  a  number  of 
(finite  dimensional  realizations  of)  p.d.e.’s.  In  practice,  the  vectors  i>i,  i  =  1  ■  ■  ■  N, 
could  be  the  n—  vector  V(r),  i,-),  when  f)  is  an  array  of  n  points  distributed  in 
space  and  V  is  some  flow  component  which  is  measured  at  the  sequence  of  times 
ti,i  =  1, • • • iV.  The  (i,j)  element  in  the  inertia  tensor  is  then  the  correlation 
function  <  V(ri,t)V(fj,t)  >.  A  major  reason  for  the  importance  of  this  basis  is 
that  if  v(t)  is  a  generic  solution  in  A  and  2,  •  •  •  is  any  orthogonal  basis  for  H, 
then  the  r.m.s.  error  of  a  k- th  order  truncation 


lim  T"1 

T— *00 


/ 


II  v(t)  -  £(»(()  •  h)h  II2  dt 

1= 0 


(13) 


is  minimal  for  all  k  >  0  precisely  "'hen  <f>j  =  tj  for  all  j  >  1. 

The  local  version  of  this  procedure  is  equally  important,  because  of  the  intermit- 
tency  and  fluctuation  of  the  organized  structures.  Like  the  filaments  of  Langmuir 
turbulence,  they  will  often  have  a  transient  character  and  averages  taken  over  all 
time  or  over  the  global  attractor  will  miss  them.  Accordingly,  we  define  the  notion 
of  a  local  SVD  as  follows.  Given  a ;  €  A,  let  Bt(x)  denote  the  ball  of  radius  e  about 
x.  Instead  of  taking  the  whole  of  A,  we  regard  the  points  in  B((x)  as  making  up 
the  mass  with  their  weights  again  given  by  u,  but  suitably  normalized.  Just  as 
in  the  global  case,  we  can  then  obtain  am  orthogonal  decomposition  of  H  and  the 
associated  moments.  We  assume  that,  as  e  — >  0,the  axes  converge  to  the  axes 

e(x):  e!(x),e2(x),e3(x),... 


and  have  the  spectrum 


a(x)  :  <Tj(x)  >  cr2(x)  >  cr3(x)  >  •  •  • 


If  the  flow  field  is  locally  dominated  by  a  d-dimensional  structure,  then  the  spec¬ 
trum  of  moments  cr(x)  has  the  property  that  the  <jj(x )  decay  exponentially  in  j  for 
j  >  a.  This  feature  is  very  important  as  it  allows  us  to  contruct  a  local  analogue  of 
the  Lyapunov  spectrum,  which  is  defined  as  the  set  of  growth  rates  of  successively 
larger  sub  volumes  of  the  tangent  space  of  a  trajectory  in  the  phase  space.  Under 
fairly  weak  assumptions  (sufficient  hyperbolicity  on  the  attractor  to  ensure  that 
any  trajectory  eventually  covers  the  attractor),  the  Lyapunov  spectrum  exists  as 
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a  global  quantity.  It  would  be  tempting  to  suggest  (hope!)  that  this  is  a  natural 
decomposition  of  the  Lyapunov  spectrum  into  a  small  set  of  a  d  order  one  exponents 
and  a  much  larger  set  D  whose  magnitudes  cluster  around  zero.  In  general,  this 
will  not  be  true,  because  the  homoclinic  excursions  are  dramatic  local  events  which, 
although  responsible  for  transport,  may  be  relatively  rare  and  therefore  a  global 
average  using  the  usual  density  of  points  invariant  measure  on  the  attractor  will 
smear  local  large  values  of  the  growth  rate  over  the  whole  trajectory.  Therefore  it 
is  unlikely  that  the  Lyapunov  spectrum  will  decompose  as  one  might  optimistically 
hope.  Neither  can  the  Lyapunov  exponent  be  defined  locally.  However  the  local 
SVD  decomposition  is  well  defined  and  does  discriminate  between  various  directions 
in  the  local  tangent  space.  For  example,  if  we  take  a  ball  on  the  attractor  at  a  point 
x  from  which  a  the  homoclinic  excursion  begins,  the  d  directions  associated  with 
the  dominant  moments  inertia  crj(x)  do  span  what  we  have  called  the  order  one 
unstable  manifold  of  Ah  at  x. 

Moreover,  we  can  adapt  these  further  in  order  to  look  for  correlations  between 
organized  structures  and  the  transport  properties  of  the  flow.  To  do  this  we  have 
to  bring  in  the  ideas  of  Broomhead  and  King  on  time-series  analysis.  Suppose  that 
besides  measuring  a  representation  of  the  state  v(t)  at  time  i,  we  also  record  a  mea¬ 
sure  n(t )  of  the  transported  quantity.  For  simplicity  of  notation  and  compatibility 
with  real  data  we  assume  that  the  time  t  is  discrete.  For  each  t  consider  the  vectors 

w(t)  =  (v(f),  n(t  -  a), . . . ,  n(i  -  6)). 


Let  A  denote  the  set  fl^x)  {u>(t)}  and  construct  the  axes  e,  and  moments  of 
inertia  cr,  for  A  as  for  the  attractor  A  in  the  global  procedure  above.  In  practice 
A  will  be  finite  so  for  v  one  takes  the  measure  on  A  which  give  the  points  equal 
weight.  For  small  i  the  axes 


fi«  —  (®i>  71— a  >  •  •  •  >  ftfc) 


will  then  give  the  dominant  structures,  each  being  an  important  spatial  structure 
e«  together  with  the  transport  time-series  n_a,...nj  correlated  with  this.  It  thus 
gives  us  the  relationship  between  spatial  structures  and  transport.  The  problem 
is  that  the  spatial  structures  which  transport  the  most  may  not  correspond  to  the 
dominant  moments  because  of  their  infrequent  occurrence  (and  hence  low  weight 
with  respect  to  v),  or  because  they  may  get  mixed  (as  linear  combinations)  with 
modes  transporting  less.  To  overcome  this  problem  we  propose  that  the  vectors  v{t ) 
should  be  weighted  by  A n(t  -j-  c)  for  some  suitable  c  and  A  >  0.  Provided  A  is  large 
enough,  the  dominant  moments  should  then  correspond  to  those  which  transport 
most.  For  example,  A  might  be  chosen  to  be  the  transport  itself  or  some  suitable 
power  thereof.  As  well  as  obtaining  the  relationship  between  spatial  modes  and 
transport  properties,  this  procedure  singles  out  those  modes  which  transport  most. 
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Finally,  we  note  that  one  can  perform  a  local  version  of  the  global  axes  of  inertia. 
This  may  be  especially  useful  when  there  are  problems  with  phase  fluctuations  or 
variations  [20]. 
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Figure  Captions 


(1)  Snapshots  of  |?/>(x)|2  at  times  before  (a)  during  (b)  and  after  (c)  the  collapse 
of  a  filament  in  the  turbulent  ensemble  sustained  by  driving  with  i/j  =  0.2. 
The  nearly  simultaneous  collapse  of  two  filaments  (d)  and  their  remnants  (e 
and  f)  shortly  after  collapse.  Note  the  outward  propagation  and  diminution 
of  the  concentric  cylindrical  shells  that  comprise  the  remnants. 

(2)  (a)  The  instantaneous  rate  of  dissipation  7(t)  =  2  7e(&)|V>(il)|2,  its  integral 

r(t)  =  —  f*  i(tf)dtf,  and  the  global  spatial  maximum  of  |^>(x)|2  ^  functions 
of  time  in  the  strong  turbulent  regime  sustained  with  =  2. 

(b)  As  in  (a)  but  for  a  longer  time.  Notice  that  strong,  sudden  dissipation 
is  exactly  correlated  with  collapse  and  that  the  loss  of  energy  (— T)  is  approxi¬ 
mately  linear  over  long  times. 


(3)  The  total  energy  dissipated  by  collapsed  filaments  ( 8nc )  divided  by  the  total 
energy  dissipated  (Sri?)  for  three  different  turbulent  ensembles  differing  in 
mean  intensity,  <  >,  as  a  function  of  the  cut-off  rate  70  used  to  define 

filament  collapse.  (See  text  for  discussion.)  Note  that  this  ratio  necessarily 
approaches  one  as  70  — ►  0. 


(4)  A  histogram  of  energies  ( 6n )  dissipated  in  the  strong  turbulent  regime  main¬ 
tained  by  driving  with  =  2.  The  cut  off  delta.  70  =  10.  Note  that  the 
distribution  shifts  towards  the  critical  value  p  =  .29  as  70  -+  0  (see  texty  for 
discussion). 

(5)  (a)  As  in  Fig.  2(a)  but  confined  to  a  small  domain  or  “box”  centered  at  a 
single  collapse  site  in  the  case  75  =  .2.  (The  time  origin  has  been  redefined 
for  convenience.) 


(b)  The  total  energy  in  the  box  ( n&ox  =  Jbox  d  — >  x|^(— >  x)\2  and  the  conser¬ 
vative  flux  through  the  boundary  of  the  box 

(flux  =  f  dt'2Im  f  dlh  •  (ip  v  ^*)) 

'  J  J  dbox 

as  functions  0^  time.  Note  the  “shoulder”  in  the  n\,ox  graph,  corresponding  to 
the  fluxing  of  1  ■»  collapse  remnant  out  of  the  box  immediately  following  the 
sudden  burnout  of  the  core  of  the  filament.  Adding  the  energy  lost  in  burnout 
(i  =  .08)  to  that  carried  away  by  the  remnant  (.08  <  t  <  .12,  i.e.,  the  shoulder) 
we  obtain  very  nearly  the  critical  value  v  —  .29. 


(6)  A  histogram  of  times  passed  between  successive  collapse  events  in  the  case 
Vf,  —  2.  Note  that  the  distribution  is  nearly  Poisson,  with  a  mean  time 
between  events  of  <  r  >=  .08.  Fluctuations  are  still  affecting  the  tail. 
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(7)  Mean  dissipation  rate  <  7  >=  J2k  27«(&)  <  \4>(k)\2  >  and  mean  time  between 

collapse  as  a  function  of  the  level  of  turbulence  (<  >  for  several  different 

cases. 

(8)  Flow  trajectories  in  the  turbulent  regime  with  r/&  =  .85  covering  about  10 
collapse  events.  The  trajectories  are  projected  onto  the  global  observables 

|4j^aa?)  to  produce  a  curve  in  3  dimensions  displayed  in  (a).  In  (b) 
we  have  projected  this  curve  onto  the  three  coordinate  planes  to  aid  visual¬ 
ization.  Those  trajectories  that  are  fluctutions  about  homoclinic  excursions 
are  the  smooth  sparse  loops  passing  through  the  larger  values  of  \ip\2max. 
The  denser,  more  jittery  part  of  the  curve  lies  in  A}j- 

(9)  Angle-averaged  correlation  function. 

(10)  The  damping  function,  7 e(k).  Note  the  rapid  turn-on  of  damping  at  k0  = 
30  =  .2 kd,  where  k&  =  | \/M  and  in  all  cases  studied  M  =  7344. 

(11)  A  contour  plot  of  the  long-time  averaged  spectral  energy  density  <  |^>(fc)|2  >. 
The  scale  is  logarithmic:  energy  densities  on  adjacent  solid  contours  differ  by 
a  factor  or  100.  Note  that  the  turbulence  is  isotropic  and  that  the  damping 
(Fig.  3)  confines  the  (time-averaged)  energy  to  long  and  intermediate  scales. 

(12)  A  schematic  drawing  of  the  attractor  and  the  homoclinic  excursion. 

(13)  Phase  plane  for  x"  —  -x  +  2x3. 

(14)  A  scenario  for  the  homoclinic  excursion  associated  with  turbulent  spots. 
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S'creamwise  Evolution  of  Naturally  Occurring  Instabil  i  ti  es 


R.  A.  Petersen 

.  Power  spectral  densities  of  streamwiss  component  u'  measured 

at  r/D  =  0.5  and  at  various  streamwise  positions  within  the 

laminar/turbuient  transition  region  of  aiii symmetric  jet  at  Li _ j  = 

16  m/s.  The  initial  instability  -frequency  is  near  1200  Hz. 

/  ■ 

Successive  parings  resul t  in  spectral  peaks  at  60S  Hz  and  300 
Hz.  The  inertial  subrange  5/3*)  emerges  near  x/D  =  0.8. 

Fourier  Decomposition 

Since  the  velocity  -field  is  statistically  homogeneous  and 
periodic  in  the  aximuthal  direction,  it  is  appropriate  to 
decompose  the  -field  into  Fourier  modes.  Measurements  at 
Cx/D,r/D3  =  C0.2,0.53,  jet  speed  =  16  m/s. 

TOR:  Sample  .time  series  from  eight  hot  wire  sensors.  Sample 

1024  paints  digitized  at  50  Hz. 


BOTTOM:  Decomposed  time  series.  Naturally  occurring 

instabilities  are  predominantly  axisymmetric. 


Phase  Portrait 

Stereoscopic  representation  phase  portrait  embedded  in  three 
dimensions.  Modes:  aQ(tj,  b-^Ct),  a2(t);  jet  speed  =  16  m/s; 
measurement  location  Cx/D,r/Dl  =  110.2,0.53;  1024  points  2  50  kHz. 

Even  though  the  time  series  exhibits  temporal  and  spatial 
coherence,  the  phase  portrait  appears  high  dimensional. 
Consequently,  an  attempt  is  made  to  develop  a  statistical  model 
from  the  phase  portrait. 
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Phase  Space  Statistics 

Based  on  the  -fact  that  the  Wavier  Stokes  equations  involve 
first  differences  in  time  and  that  the  (incompressible)  pressure 
is  determined  by  the  instantaneous  velocity  field,  we  model  the 
system  as  a  noise  driven  Harkov  process.  That  is.  the 
expectation  of  the  future  state  of  the  system  depends  only  on  the 
present . 

Tne  conditional  mean  represents  the  temporally  coherent  part 
of  the  decomposed  signal  and  is  related  to  the  temporal  cross 
correlation  function.  The  "fine  scale-'1'  turbulence  is  modelled 
by  a  random  vector  of  zero  mean  and  unity  variance.  The 
conditional  covariance  can  be  decomposed  locally  into  a  set  of 
eigenvectors  and  eigenvalues.  If  the  conditional,  joint 
statistics  are  uncorrelated  one  would  expect  the  conditional 
eigenvectors  to  parallel  the  basis  vectors  of  the  original 
decomposi ti on . 

Measured  Conditional  Accelerations 

Joint  statistics  are  based  on  a0(t),  horizontal  axis,  and  on 
b  (t> ,  vertical  axis.  Measurements  span  the  region  from  laminar 
flow  to  fully  developed  (based  on  emergence  of  inertial  subrange) 
turbulence.  At  x/D  =  0.2  the  conditional  convariance 
eigenvectors  are  correlated  by  the  instability!!  by  x/D  =  0.8  the 
conditional  means  are  small  compared  to  the  conditional 
covariance  eigenvalues  (root  mean  square)  and  the-  conditional , 
joint  statistics  appear  to  be  uncorrel ated . 

In  each  case  the  sampling  rate  was  50  kHz,  and  the  joint 
statistics  are  based  on  1,024,000  data  points.  The  data  was 
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digitised  in  1024  point  records  and  de-meaned  before  h&i  ng 
decomposed  into  Fourier  modes.  The  display  is  modulated  by 
contours  of  joint  probability  density:  contours  between  50—70'/. 
and  between  80-1007,  of  the  peak  value  were  selected.  The  joint 
probability  distributions  are  evidently  skewed  and  the  display  is 
approximately  4  ;<  4  standard  deviations.  Approximate  because  the 
S'  aline)  is  based  on  I, hr-  first  17  of  the  database. 


